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Abstract 

The question of B.H. Neumann, which dates back to the 1950s, asks 
if there exists an outer billiards system with an unbounded orbit. We 
prove that outer billiards for the Penrose kite, the convex quadrilateral 
from the Penrose tiling, has an unbounded orbit. We also analyze some 
finer properties of the orbit structure, and in particular produce an 
uncountable family of unbounded orbits. Our methods relate outer 
billiards on the Penrose kite to polygon exchange maps, arithmetic 
dynamics, and self-similar tilings. 

1 Introduction 

1.1 History of the Problem 

Outer billiards is a basic dymamical system which serves as a toy model for 
celestial mechanics. See Sergei Tabachnikov's book [T], and also the survey 
[DT], for an exposition of the subject and many references. 

To define an outer billiards system, one starts with a bounded convex set 
S C R 2 and considers a point x E R 2 — S. One defines X\ to be the point 
such that the segment xqX\ is tangent to S at its midpoint and S lies to 
the right of the ray xqx{. (See Figure 1.1 below.) The point x\ is not well- 
defined if XqX\ is tangent to S along a segment. This will inevitably happen 
sometimes when S is a polygon. Nonetheless, the outer billiards construction 
is almost everywhere well defined. The iteration xq — > X\ — > a^--- is called the 
forwards outer billiards orbit of Xq. The backwards orbit is denned similarly. 



* This research is supported by N.S.F. Grant DMS-0604426 
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Figure 1.1: Outer Billiards 

B.H. Neumann introduced outer billiards during some lectures for pop- 
ular audiences given in the 1950s. See, e.g. [N]. J. Moser popularized the 
construction in the 1970s. Moser [M, p. 11] attributes the following ques- 
tion to Neumann circa 1960, though it is sometimes called Moser's Question. 

Question: Is there an outer billiards system with an unbounded orbit? 

There have been several results related to this question. 

• Moser [M] sketches a proof, inspired by K.A.M. theory, that outer 
billiards on S has all bounded orbits provided that OS is at least C 6 
smooth and positively curved. R. Douady gives a complete proof in his 
thesis, [D]. See [B] for related work. 

• In [VS], [Ko], and (later, but with different methods) [GS], it is proved 
that outer billiards on a quasirational polygon has all orbits bounded. 
This class of polygons includes polygons with rational vertices and also 
regular polygons. In the rational case, all defined orbits are periodic. 

• Tabachnikov analyzes the outer billiards system for the regular pen- 
tagon and shows that there are some non-periodic (but bounded) or- 
bits. See [T, p 158] and the references there. 

• Genin [G] shows that all orbits are bounded for the outer billiards 
systems associated to trapezoids. He also makes a brief numerical study 
of a particular irrational kite based on the square root of 2, observes 
possibly unbounded orbits, and indeed conjectures that this is the case. 

1 My information on this comes from 1999 email correspondence between Bernhard 
Neumann and Keith Burns, and also from Bcrhnard's son Walter. 
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1.2 The Main Result 



The main goal of this paper is to prove: 

Theorem 1.1 Outer billiards on the Penrose kite has an unbounded orbit. 
In fact, both the forwards and backwards orbits of the point p shown in Figure 
1.2 are entirely defined and unbounded. 




Figure 1.2: The Penrose Kite 



The Penrose kite is the convex quadrilateral that appears in the Penrose 
tiling. (In general, a kite is a convex quadrilateral with bilateral symmetry.) 
Figure 1.2 shows a classic construction of the Penrose kite-the shaded figure- 
based on a regular pentagon. The additional lines show how the point p is 
constructed. The significance of the points v and w will be explained in 
Theorem 11.21 below. In §2.1 we will give more traditional coordinates for the 
objects of interest to us. 



3 



1.3 Outline of the Proof 



We will give a rigorous computer-assisted proof of Theorem 11.11 We think 
that our proof gets close to the conceptual core of what is going on, but we 
need a lot of computation to make it go. Here are the main ideas. Let T 
denote the square of the outer billiards map. 

1. We replace the Penrose kite by an affinely equivalent kite S whose 
vertices lie in the ring Z [(/>], where <fi is the golden ratio. We show that 
the outer billiards map (and in particular T) is entirely defined on a set 
C(±), consisting of certain points whose first coordinates are positive 
and lie in \Z[(f)] and whose second coordinates are ±1. See Figure 2.1. 
Let Tfl : C(±) — > C(±) be the first return map of T. 

2. Forgetting about the ^/-coordinate (which is always ±1) we identify 
C(±) with the set of vertices in Z 2 fl H, where H C R 2 is a certain 
halfplane. We then encode the dynamics of Y# by a graph T C H. 
The vertices of T lie in Z 2 fl H and the edges are short lattice vectors. 
We call T the arithmetic graph. See Figures 2.2 and 2.3 for pictures. 

3. The structure of V is controlled by a certain partition V of the square 
torus T 2 into 26 convex polygons, and an associated dynamical system 
akin to a polygon exchange map. See Figure 2.5. This structural result, 
which we call the Arithmetic Graph Lemma, is the central technical 
result in the paper. 

4. The dynamical process described in Item 3 is compatible with a kind of 
multi-valued contraction, defined on various simply connected subsets 
of T 2 . The graph of this correspondence is an irrationally embedded 
complex line sitting inside T x T . Rick Kenyon tells me that this 
picture is essentially the same as what one sees in the cut and project 
method from the theory of quasi-crystals. 

5. The compatibility discussed in Item 4 forces V to behave like an ape- 
riodic tiling with an inflation rule. Compare [Ke]. That is, when we 
dilate V about the origin by 3 we find that the dilated image is very 
closely shadowed by the original image. This structure forces the con- 
nected component Tq of V containing the origin to travel unboundedly 
far away from OH, and this translates into the unboundedness of the 
orbit of our point p shown in Figure 1.2. 



4 



1.4 Orbit Structure 



Our proof of Theorem 1 1 . 1 1 gives us information about the fine orbit structure 
of outer billiards on the Penrose kite. We will see (Lemma 12. 5j) that the 
Euclidean norm of the nth point of + (p) is a proper function of n. Here 
+ (p) is the forward orbit of p. (A function / : N — > N is proper if n] 
is bounded for all n G N.) In contrast, the backwards orbit 0_(p) returns 
densely to a certain Cantor set. 

Given a horizontal interval / C R 2 and some A G (0, 1) let C(I, X) C I 
denote the (Cantor set) limit set of the semigroup generated by the two 
contractions of strength A that map I into itself and fix an endpoint of /. 
We say that a gap point of C = C(A, I) is an endpoint of any component of 
I — C . We let C* denote C minus its countably many gap points. 

Let 70 (respectively 71) be the similarity of strength —A (reversing orien- 
tation on J) that maps C(A, /) to its own right (respectively left) half. Let Z 2 
denote the 2-adic integers. There is a unique homeomorphism 9 2 : Z 2 — > C 
that conjugates jj to the map x —* 2x + j. Let p, v, w be as in Figure 1.2. 
We have set things up so that #2(0) = p when / = [v, w] and A = (f)~ 3 . 

Theorem 1.2 Let C = C([v , w], <p~ 3 ), where v and w are as in Figure 1.2 
and (j) is the golden ratio. The forwards and backwards orbits of every point 
of C* are defined and unbounded. 

• If x G C* — #2(0) then the T -forwards orbit of x first returns to C* as 
X-, such that 2 ~ 1 {x-) = — 1- The number of iterates between x 
and x_ , as well as the maximum Euclidean norm of such an iterate, is 
a proper function ofl/\\x — ^2(0) || . 

• //i6C*-^(-l) then the T -backwards orbit of x first returns to C* 
as x + , such that ^ 2 ~ 1 ( a; +) = + 1- The number of iterates between 
x and j CIS well as the maximum Euclidean norm of such an iterate, 
is a proper function of l/\\x — 9 2 {— 1)|| . 

In the two items, we mean to consider just the even iterates in the orbits. 
Again, T is the square of the outer billiards map. Theorem 11.21 immediately 
implies that there are uncountably many unbounded orbits of the outer bil- 
liards map on the Penrose kite with the following wild behavior: If V is any 
neighborhood of the vertex v and U is any neighborhood of 00, then both the 
forwards and backwards orbits oscillate between U and V infinitely often. 
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1.5 Discussion 



Here are some remarks on the origins of Theorem 11.11 Let 5" be the Penrose 
kite. It appears that R 2 - S' = AUBUC where 

• A is a dynamically invariant union of finite sided polygons. All the 
orbits in a polygon are periodic with the same combinatorial type. 

• B is a countable union of line segments consisting of the points on which 
the outer billiard map is undefined. B is the so-called discontinuity set. 

• C is a fractal set consisting entirely of unbounded orbits. 

Tabachnikov develops a similar picture for the regular pentagon, except 
that his region C consists of non-periodic but bounded orbits. The general 
regular iV-gon seems to have a similar kind of structure, though nobody has 
yet made a detailed study. Tabachnikov has a beautiful picture of the case 
N = 7 on the cover of his book [T]. Inspired by Tabachnikov's picture, I 
designed a computer program, Billiard King, d which draws these special 
sets for kites and explores their structure. 

When the kite has rational vertices, C is empty and A is a locally finite 
tiling. I searched B through the parameter space of rational kites for examples 
where A featured small tiles having widely ranging orbits. Eventually I 
considered kites having vertices with Fibonacci number coordinates, and this 
led to the Penrose kite. 

I hope to find purely conceptual proofs for the results in this paper, but 
my first goal is just to establish them as true statements. Many people do 
not like computer-aided proofs, so I would like to say several things in favor 
of the proof I have given. First of all, it is the best I could do and there 
are no other proofs for results like these. Second, the proof is not just a 
calculation— there are plenty of concepts in it. Finally, I was able to check 
essentially all the details in the proof using Billiard King, so this proof has a 
lot of computational safeguards that one might not see in a traditional proof. 

I tried as hard as possible to write a proof that is independent from 
Billiard King, but still the reader of this paper would get a much greater 
appreciation for what is going on by learning to use Billiard King. A few 
minutes playing with Billiard King is worth hours of reading the paper. 

2 Onc can download this Java program from my website: www.math.brown.edu/~res 
3 Had I known earlier about Genin's numerical study, mentioned above, I perhaps would 
have saved time. 
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1.6 Organization of the Paper 

Part 1: (§2-5) In this part we prove all our main results modulo the Arith- 
metic Graph Lemma. (See Item 3 of §1.3.) In §2 we reduce Theorem 11.11 to 
the statement that To, the component of the arithmetic graph of interest to 
us, admits what we call an inflation structure. Essentially this means that 
certain smallish subsets (which we call genes) of T are, when dilated by <p 3 , 
closely shadowed by other subsets of r . In §3 we tie the existence of an 
inflation structure to the dynamical behavior of the torus partition guaran- 
teed by the Arithmetic Graph Lemma, and thereby reduce the problem to a 
finite calculation with integer arithmetic. In §4 we explain how we perform 
the calculation. In §5, we prove Theorem 11.21 by taking a close look at the 
specifics of our inflation structure. 

Part 2: (§6) In this part, we prove the Arithmetic Graph Lemma. To 
this end, we factor the return map T# (from Item 2 of §1.3) as the product 
of 8 simpler maps, which we call strip maps. (Actually, the decomposition 
involves a 9th map as well. See Equation [321 for a precise statement.) We 
then show how each strip map actually is compatible with a certain embed- 
ding of R 2 - S dense subspace of the 4-torus. Each strip map extends 
to act as a piecewise affine transformation of the 4-torus. We then embed 
the dynamics of as a 2-dimensional geodesic slice of this 4-dimensional 
dynamical system. The Arithmetic Graph Lemma is a consequence of this 
structure and some additional integer arithmetic calculations. 

Part 3: (§7-8) In §7 we describe our calculations using pesudo-code, in 
enough detail that the interested reader should be able to reproduce them. 
In §8 we include the data used in our calculations. All the calculations are 
implemented in Billiard King, a Java program. Billiard King is a massive pro- 
gram, but we have placed the smallish number of routines actually relevant 
to the proof in separate files so that they are easier to survey. 

1.7 Acknowledgements 

I thank Jeff Brock, Keith Burns, Peter Doyle, David Dumas, Eugene Gutkin, 
Pat Hooper, Richard Kent, Rick Kenyon, Curt McMullen, and Steve Miller 
for helpful and interesting discussions related to this work. I especially thank 
Sergei Tabachnikov for all his help and encouragement. 
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2 The Proof in Broad Strokes 



2.1 An Affinely Equivalent Shape 

The Penrose kite is affinely equivalent to the quadrilateral with coordinates 
(0,1); (-1,0); (0,-1); (A,0); A = 0~ 3 (1) 

where <fi is the golden ratio. The affine equivalence maps the point p in Figure 

1.2 to the point 

(<r 2 ,-i)- (2) 

To clarify the logic of our argument, it is useful to consider a general 
value of A e (0, 1) for the moment. We define T : Z 2 — > R via the formula 

T(x,y) = 2Ax + 2y + ^-. (3) 

The point in Equation [2] is (T(0, 0), —1) when A = 0~ 3 . 
Lemma 2.1 TTie outer billiards orbit of any point in the set 

T(Z 2 ) x {±1} 

is entirely defined, both forwards and backwards. In particular, the entire 
orbit of the point in Equation [H is defined. 

Proof: Let L denote the family of horizontal lines in R 2 whose y-coordinates 
are odd integers. The outer billiards map preserves L. The only points where 
the first iterate of the outer billiards map is undefined are points of the form 
/ n e, where I is a line of L and e is a line extending an edge of our kite. 
One can check easily, given Equation [TJ that all such points have first coor- 
dinates of the form m + An, where m, n G Z. On the other hand, no point 
of T(Z 2 ) x {±1} has this form, and no iterate of such a point has this form. 4k 

Let T denote the square of the outer billiards map. Let C = T(Z 2 )x Z 0( jd, 
where Z odd is the set of odd integers. The vector T(x) —x is always twice the 
difference between two of the vertices of our shape S. Given Equation II. 1\ 
we see that the first coordinate of T(x) — x always has the form 2m + 2nA 
where m and n are integers. The second coordinate is always an even integer. 
Hence T is entirely defined on C and preserves this set. 
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2.2 The Arithmetic Graph 

Let 

C(±) = (T(Z 2 )n J R + )x{±l}, (4) 

where T is as in Equation [3l and A = <ft~ 3 . (Henceforth we take A = <f)~ 3 .) 
The set C(±) is dense in the union of two rays starting at (0, ±1) and parallel 
to (1,0). The point in Equation [2] belongs to C(-). 




Figure 2.1: The objects S and C(±). 
T does not preserve C(±), but we have0 the first return map: 

T« : C(±) - C(±). (5) 

Our idea is to encode the structure of T R graphically. We join Xi,x 2 G Z 2 
by a segment iff T(xj) > for j = 1, 2, and there are choices e±, e 2 G {—1, 1} 
such that 

T R (T(x 1 ),e 1 ) = (T(x 2 ),e 2 ). (6) 

There is a subtlety in our definition, caused by the asymmetric roles played 
by X\ and x 2 in the construction. One might worry that this asymmetry could 
lead to an inconsistency, whereby we are told to join Xi,x 2 by an edge, but 
then told not to join x 2 , x\ by an edge! However, we observe that reflection 
in the X-axis interchanges C(+) and C(— ) and conjugates the outer billiards 
map to its inverse. Therefore 

T R (T(xi), ei) = (T(x 2 ), e 2 ) T R (T(x 2 ), -ei) = (T(xi), -e 2 ). 

Hence, our definition is consistent. 

We let r C H denote the graph resulting from our construction. Here 
H = T _1 (i2 + ) C R 2 is an open halfplane containing T. We call T the 
arithmetic graph. 

Figure 2.2 shows a small portion of T. The dot is (0, 0). The component 
To containing (0, 0) is drawn in black. The grid indicates Z 2 . The black line 
running along the bottom is dH. 

4 A "runway" like the one in Figure 2.1, as well as a map somewhat like T#, is also 
considered in [GS]. 
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Figure 2.2: Part of the arithmetic graph 




Figure 2.3: More of the component IV 



Figure 2.3 shows more of To, again drawn in black. The origin is denoted 
by a little vertical line segment touching OH at its top endpoint. We have 
erased the other components and deleted the grid, to get a clearer picture. 
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(The reader should use Billiard King to see really great pictures of T.) The 
curve T , and indeed all of T, behaves like the self-similar inflationary tilings 
studied in [Ke] . One might call T a "large-scale fractal" . 

Here is what we mean by inflationary: In Figure 2.3, there is also a grey 
curve running alongside IV This grey curve is the dilated image 3 IV In 
other words, we dilate T by </> 3 (about the origin), color it grey, and superim- 
pose it on the original picture. From what we can see of the picture, it looks 
like To and 3 Fo closely follow each other. One might say that Tq is quasi- 
invariant under a dilation. If this is true— and we will prove it below— then 
both ends of To must exit every tubular neighborhood of dH. (Below we 
formulate this principle precisely.) Theorem 11.11 follows immediately from 
the fact that both ends of T exit every tubular neighborhood of dH. 

2.3 The Arithmetic Graph Lemma 

Let T 2 = (R/ Z) 2 be the square torus. Given p E R 2 let [p] denote the 
projection to T 2 . We define \1/ : Z 2 — > T 2 with the equation 



<t>{x,y) 



' T{x,y) T(x,y) 
v 20 ' 2 



*x + < 



v| 4 



(7) 



Here T is the map from Equation [3j The second equation is a direct calcu- 
lation, which we omit. ^(Z 2 ) is dense in T 2 . 

Given v G Z 2 D H we define the local type of v to be the translation 
equivalence class of the union of edges of T emanating from v. It turns out 
that there are 16 types (including the type where no edges emanate from v) 
but we prefer to label these 16 types by the integers 1, ...,26 according to 
Figure 2.4. Most of these types can be seen in Figures 2.2 and 2.3. 



5 7 



6 8 



10 



12 14 


15 17 


16 18 


















/ 













20 21-23 24-26 



Figure 2.4: The local types 
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Figure 2.5 shows a partition of T 2 into 26 open convex polygons Pi, P26- 
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/7 19/ 
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7 / 




/ 21 


\ / 12 \ 



Figure 2.5: The Torus Partition 



We have changed the aspect ratio so as to get a nicer picture. The true 
geometric picture is contained in the unit square, with sides identified as 
usual. The origin (0, 0) is the bottom left corner. The only polygon that 
"wraps" around is P 2 6- In §6 we will prove 

Lemma 2.2 (Arithmetic Graph) ^ maps each x G Z 2 n H into some 
open polygon ofV, and x has local type k if and only if^(x) G P&. 

The next lemma is not needed for our theorems, but it is a nice fact to 
know and it illustrates the power of the Arithmetic Graph Lemma. 

Lemma 2.3 The arithmetic graph is an embedded union of polygons and 
polygonal arcs. 

Proof: (Sketch.) First of all, we check the list of types that appear in the 
Arithmetic Graph Lemma and we observe that every vertex in the graph 
has valence 2. Thus, the only kind of crossing that can occur is a transver- 
sal crossing, where the cross point is not in Z 2 . In particular, there would 
have to be a point (x,y) of type a and a point (x + l,y) of type b, where 
a G {4,6,8} and b G {2,4,9, 12, 14}. But v = V(x + l,y)-if?(x,y) is the vec- 
tor (0" 3 ,0 -4 ). This vector is the one that points from the lower left corner 
of P 3 to the upper right corner or P 3 in Figure 2.5, when it is anchored at 
the point (0, 0). Looking carefully at Figure 2.5 we see that no position of v 
has this property. 4 
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2.4 Inflation Structures 



As we mentioned above, the basic idea in our proof that T rises unboundedly 
far away from dH is to show that r is quasi-invariant under the dilation 

$(x,y) = (0 3 x,0 3 y). (8) 

In this section we formalize this idea. 

Say that a gene is a (combinatorial) length 6 connected component of r . 
It turns out that there are 75 genes up to translation equivalence. Each gene 
A has a core B, consisting of the central path of length 2. The left hand 
side of Figure 2.6 shows a gene and its core. We say that p' E Z 2 shadows 
p £ R 2 if \\p — p'\\ < 4. (Here 4 is a convenient cutoff.) Let B be a gene core, 
as above. We say that a finite polygonal path A' C T shadows Q(B) if each 
endpoint of A' shadows a corresponding endpoint of Q(B). Figure 2.6 shows 
a combinatorially accurate example. 




L^J-J 1 L_ 

Figure 2.6 genes, gene cores, and shadowing 



Let Q denote the set of all genes. Let X denote the set of all finite polygo- 
nal paths of T. We say that an inflation structure is an assignment x, to each 
gene A e Q a path A' = x(A) in X which shadows Q(B). We emphasize that 
x{A) only shadows the dilated core $(-£>), but perhaps depends on all of A 
for its definition. We also emphasize that the translation type of x(A) only 
depends on the gene type of A. Billiard King computes part of an inflation 
structure and allows the interested user to interact with it. 

We say that a inflation structure is coherent if 
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• If A is the gene whose center vertex is (0, 0), then A C x{A)- We call 
A the zeroth gene. 

• If A\ and A 2 are two consecutive genes, then x(Al) U x(A 2 ) is a sin- 
gle polygonal path. In other words, an d xi^) splice together 
seamlessly, just as Ay and A 2 splice together. 

Lemma 2.4 7/T admits a coherent inflation structure then either direction 
of Tq rises unboundedly far away from dH, and hence Theorem \l.l\ is true. 

Proof: We see explicitly by direct computation that either direction of r 
contains a polygonal arc which starts at (0, 0) and rises 100 units away from 
dH. Let Ti be one of these two directions. Say that we have already shown 
that Ti rises up d > 100 units away from dH. Let Aq, Ay, A 2 , ... be the con- 
secutive genes of r 1; rising up to this height, chosen so that A is the zeroth 
gene. Then A C x(A ) and the consecutive paths x(^-o), xi^),... 
fit together. These paths therefore trace out Ti and rise up at least (say) 
4> 3 d — 100, a quantity larger than 2d. Thus, we have now shown that Ti rises 
up 2d units from dH in the same direction. Iterating, we get our result. A 

Lemma 2.5 If T admits a coherent inflation structure then the distance 
from the Nth point of + (p) to the origin is a proper function of N. 

Proof: Let Tq be the component of Tq — (0, 0) whose first (say) 100 iterates 
are contained in the positive quadrant E? + . (See Figure 2.3.) An argument 
just like the one given for Lemma 12.41 shows that Tq C B? + . As we trace 
out Tq, a non-periodic curve with integer vertices, we can only return to 
any given compact neighborhood of (0, 0) finitely many times. Hence, the 
distance from the iVth vertex of Tq to (0, 0) is a proper function of N. Since 
C B? + and dH has negative slope, we now see that the distance from the 
iVth vertex of Tq to dH is a proper function of N. Hence the Euclidean norm 
of the Nth. point of + (p) DC(±) is a proper function of N. But now observe 
that the norms of the points in 0+(p) between two consecutive points of 
O+ip) H C(±) are, up to a uniformly bounded factor, the same as the norms 
as the two points of + (p) flC(±). (See Figure 6.1. and the Pinwheel Lemma 
of §6.2.) ^ 
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3 Manufacturing an Inflation Structure 



In this chapter we assume the Arithmetic Graph Lemma and use it to delve 
more deeply into the connection between the arithmetic graph T and the 
torus partition V . 

3.1 Pointed Strands and Dynamical Polygons 

Say that a •pointed strand of T is a pair (X, x) , where X is a finite polygonal 
arc of T having length at least 2, and X IS cL vertex of X. Say that two 
pointed strands (Xi,X\) and (X 2 ,x 2 ) are equivalent if there is a translation 
of Z 2 carrying one to the other, and ^(xi), ^(x 2 ) belong to the same polygon 
Pj of V . Say that a pointed strand type of length n is an equivalence class 
of pointed strands, where the strands have n segments. Let S CT denote those 
x e Z 2 HH such that a = [(X, x)). 

Lemma 3.1 (Dynamical Polygon) There exists an open convex polygon 
P a C T 2 such that a E S CT if and only if^(a) E P a . 

Proof: We first consider the effect of keeping the strand but changing the 
basepoint. Suppose o\ = [(X, xi)] and a 2 = [(A, a^)], where x% and x 2 are 
adjacent vertices of X. We have x 2 = x± + (ei,e 2 ) where ej E { — 1,0,1}. 
Equation [7] gives us ^/(x 2 ) = +C, where C only depends on e\ and e 2 . 

If we have managed to prove this lemma for <j\ then we define P a2 = P ai + C 
and the lemma follows for er 2 . We call this the vertex slide argument. 

By the Arithmetic Graph Lemma, this lemma is true for the local types. 
The vertex slide argument now establishes this lemma for all pointed strand 
types of length 2. Now suppose that a = [(X,x)] has length at least 3. By 
the vertex slide argument, we can assume that x is an interior vertex of X. 
Hence X = X\ U X 2 where x E Xj and Xj is a shorter strand than X. By 
induction there are open convex polygons P ai and P U2 such that a E iff 
\I/(a) E P ay Let P a = P ai fl P a , 2 . Evidently P a is an open convex polygon. 
If a E S CT then a E S CT1 fl S CT2 and by induction \I/(a) E P ai fl P a2 = P a . 
Conversely suppose that ^(a) E P a . Then \l/(a) E P a and hence there is a 
strand Aj such that aj = [(Aj, a)]. Letting A = A\ U A 2 we have a = [(A, a)]. 
Hence a E T, a . This completes the induction step. 4k 

We call P a a dynamical polygon because, as we explain in §4, it can be 
produced by a dynamical process. 
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3.2 Inflation Maps 



In this section and the next we build some machinery for describing the 
interaction between the dilation map 

$(x,y) = (0 3 x,0 3 y) 

and the map $ : Z 2 — > T 2 given in Equation [7J That is, 



- 4 x + 0" 1 2/,r 3 a:) + i(^ 3 ,r 2 ) 



Our goal is to define, in a canonical way, a kind of integral approximation to 
$. Such approximations are vital to the creation of an inflation structure. 

We say that a small polygon is an open convex polygon of T 2 that is 
contained inside a square of side length 1/2. The polygons of interest to us 
will all be small. Let P be a small polygon. We say that a map 7 : P — > T 2 
is a special similarity if, relative to local Euclidean coordinates on P and 
7(P) the map 7 is just multiplication by — <p~ 3 . Put another way, 7 is a 
holomorphic map on P whose complex derivative is constantly equal to — 0~ 3 . 

Define 

E(p) = r 1 (P)nz 2 . (9) 

If P — P a then £(P) fl H = E CT , the set defined in the last section. We say 
that a map /3 : £(P) — > Z 2 is an inflation map if 

*G3(a)) = T(*(a))l Va G £(P), (10) 

for some special similarity 7 : P — > T 2 . In the next section we prove the 
existence of such maps. Here we investigate some of their abstract properties. 

Lemma 3.2 Suppose that [3 is an inflation map defined relative to P. Let (3 
be the new map defined by the equation f3(a) = (3{a) + b for some b 6 Z 2 . 
Then (3 is also an inflation map defined relative to P. 



Proof: We compute that 

*03(a)) = *G8(a)) + C = 7 (*(a)) + C = 7 (*(a)). 

Here C is some constant and 7 is some other special similarity which differs 
from 7 only by a translation. Our equation makes sense because T 2 = P 2 /A 
is canonically an abelian group. 4fc 
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Lemma 3.3 Suppose that (3 is an inflation map defined relative to P. Given 
a G Z 2 let P = P + *(a ). Define (3 : S(P)^-> Z as /3(a) = (3{a - a ). 
Then (3 is an inflation map defined relative to P. 

Proof: If a G E(P) then *(o) G P. But then *(a-a ) = *(a)-*(a ) G P. 
Thus P(a — a ) is always defined. We compute 

*09(o)) = *(/3(a - a )) = (/3(a)) -C = 7 (*(a)) - C = 7 (*(a)). 

Here C is some constant. The rest of the proof is as in the previous lemma. 4ft 

Lemma 3.4 Let f3\ and /3 2 oe t wo inflation maps defined relative to the same 
polygon P. If f3\(a) = /3 2 (a) for some a G £(P) then (3\ — (3 2 - 

Proof: Let jj be such that jj(^(a)) = ty(/3j(a)). Note that 71 and 72 must 
differ by a translation. Thus, these two maps are equal if they agree at any 
point. But 71(^(0)) = *(/3i(a)) = *(/3 2 (a)) = 72 (\I>(a)). Hence 71 = 72. 
Now we know that \I/(/3i(a)) = #(/3 2 (a)) for all a G S(P). To finish our 
proof, we note that \1/ is pretty obviously injective on Z 2 . A 

Now we lay the groundwork for explaining how (3 is an integral approxi- 
mation to <£>. Given an inflation map f3 defined relative to P and a G S(P), 
let 

Vp(a) = $(a) -/5(a). (11) 

Here vp measures the discrepancy between (3 and the dilation $. Say that (3 
is K-pseudo-Lipschitz if 

IMai) - ^(a 2 )|| < K\\V( ai ) - *(o2)||. (12) 

Lemma 3.5 Suppose that the inflation map (3, defined relative to P, is K- 
pseudo-Lipschitz. Then any inflation map (3 defined relative to P is K- 
pseudo-Lipschitz. 

Proof: In light of Lemma [3.41 we must have, for all a G S(P), the identity. 
(3{a) = (3{a) +a for some a G Z 2 . The result is obvious from this identity. 4k 
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3.3 The Shadow Lemma 

Recall that convex polygon P C T 2 is smaW if P is contained inside an open 
square of radius 1/2 in T 2 . 

Lemma 3.6 (Shadow) Let P be a small polygon. There exists an 4-pseudo- 
Lipschitz inflation map (5 defined relative to P. 

It suffices to prove the Shadow Lemma for the map 



which differs from ^ by a translation of T 2 . That is, we will produce a special 
similarity 7 such that V (P(a)) = 7(^0(0)) for all a G £(P). 

Given that P is small, there are intervals ii,/ 2 C R, each having length 
less than 1/2 with the following property: For a = (x, y) G A, we have unique 
integers m and n such that 



We write /i ~ // if /i — // G -Z. Since 3 ~ 3 and x E Z we have 



Since 3 ~ 20 1 and —20 4 ~ 30 3 and x,y E Z we have 

3 y ~ 20 _1 y = 2(0" 4 a;+0- 1 i/)-20- 4 a; ~ 2(m+ei)+30~ 3 a; ~ 2ei+3e 2 . (15) 
Accordingly we define 

P(x,y) = (x,y); x = 3 x - e 2 ; y = 3 y - 2ei - 3e 2 . (16) 
By construction (x, y) G Z 2 . 
Lemma 3.7 For any a±, a 2 G X(P) we /iawe 



4 x + 1 y = m + ei, 



3 x = n + e 2 ; 



(13) 




(14) 



u(ai) - u(a 2 )|| < 4||* (ai) - *o(a 2 )||- 
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Proof: Let (eij,e 2 j) be constants associated to the point dj = (xj,yj) as in 
Equation [T3"l Let 8j = — e 2 j for j = 1, 2. First, we have 

^o(ai) - #o(a 2 ) = (6 1 ,6 2 ). 

Second, we have 

v(aj) = $(%) - f3{aj) = (0 3 x - x, (j) 3 y - y) = (e 2j , 2e 2j + 3e 2j ). 

Hence 

v(a 1 ) - v(a 2 ) = (8 2 , 28 x + 38 2 ). 

The matrix associated to the linear map (8i,8 2 ) — > (8 2 ,28i + 38 2 ) has L 2 - 
norm equal to a/14 < 4 and hence the map itself is 4-Lipschitz. 4jk 

We compute 

0~ 3 x = 0~ 3 (0 3 x - e 2 ) 
Lemma 3.8 <f)~ A x + ~ —<j)~ 3 6i 

Proof: We have 

0" 4 x + ct>- l y = <p~ l x + <p 2 y - 0~ 4 e 2 - 2</r 1 e 1 - 3(1)-%. (18) 

But 

-1 2G + 2 y ~ _1 a; + <\T x y = - </r 4 )x + (<f)~ 4 x + ^y) = 

(2<f)~ 3 )x + (0" 4 x + ~ ei + 2e 2 . (19) 

Our last equation comes from Equations [13] and [141 Plugging Equation [19] 
into Equation [18] and grouping terms, we have 

0" 4 x + <j)- l y ~ (-20- 1 + l)ei + (-0- 4 - 30 1 + 2)e 2 = -0~ 3 ei. 

This completes the proof <jt 

Equation [T7I and Lemma I3TB1 together show that ty (x,y) = —(f)~ 3 ^o(x,y) 
in local isometric coordinates. This is what we wanted to prove. This com- 
pletes the proof of the Shadow Lemma. 



= x - (p- 3 e 2 ~ -<T 3 e 2 . (17) 
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3.4 Inflation Maps and Dynamical Polygons 

Now we put together the material from the previous two sections. The result 
we prove here, an immediate consequence of what we have already done, 
provides the key step in constructing a coherent inflation structure from a 
finite amount of computational information. 

Let <7i and a 2 be two pointed strand types. Let P ai and P a , 2 be the 
associated dynamical polygons. We write 

<J\ — >/3 cr 2 (20) 

if there is a inflation map (5, defined relative to P ai and having the property 
that 

l{P n ) C P a2 . (21) 
Here 7 is the special similarity associated to (5. 

Lemma 3.9 Suppose 07 — >@ a 2 . Suppose also that a\,a 2 G Z 2 n H are 
such that a 2 = /3(ai) and o~i = [(Ai,ai)] for some polygonal arc A\. Then 
there exists a polygonal arc A 2 C T such that a 2 = [(^2,02)]. 

Proof: We have a\ G C S(P (Tl ) by the Dynamical Polygon Lemma. 
By definition ^(02) = 7(^(01)) C P CT2 . By the Dynamical Polygon Lemma 
we have a 2 G S^. This means that there is an arc A 2 of Y such that 
o 2 = [(A 2 ,a 2 )]. 4 



3.5 The Inflation Generator 

Now we turn our attention to the genes we discussed at the end of §2. Recall 
that a gene is a polygonal arc of length 6 contained in IV Here To is the com- 
ponent of T that contains (0, 0). Each gene A gives rise to a pointed strand 
(A, a) where a is the central vertex of A. Each gene Aj therefore gives rise 
to the dynamical polygon Pj associated to the pointed strand type [(Aj, aj)}. 
It turns out that there are 75 such dynamical polygons, corresponding to 
75 distinct combinatorial types of gene. The corresponding polygons are 
the dark-shaded polygons in Figure 3.1. (We will explain the significance of 
the 24 light-shaded polygons in the next chapter.) One can see easily from 
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the picture, or from inspecting the list in §8.11 that any two points in the 
same dynamical polygon Pj are within 1/4 of each other. In particular, these 
polygons are all small. 




Figure 3.1: dynamical polygons corresponding to genes 



We say that the gene A (with central vertex a and gene core B) is nicely 
shadowed by a polygonal arc A' if there is a vertex a' of A 1 such that the 
following is true: 

• There exists an inflation map f3 with the properties that (3(a) = a' and 
P,a)]^P>')]. 

• Each endpoint of is within 3 units of an endpoint of A'. 

We say that an inflation generator is a list G of 75 pairs of the form 
{(Ai, A'j)} such that each gene of T is translation equivalent to one gene Ai 
on the list, and A-CT nicely shadows A4. 

Let A be some gene. There is some gene Aj on our list such that A 
and Aj have the same type. Let a be the central vertex of A and let aj be 
the central vertex of Aj. Let /3j, A'j, and a'j be the objects associated to 
Aj, as above. From Lemma [3.91 there is a pointed strand (A', a') such that 
[(A'j,a'j)] = [(A', a')} and a' = fy(a). We let X (A) = A'. 

Lemma 3.10 x ^ s an inflation structure. 

Proof: For an arbitrary gene A, with gene core B, we need to show that each 
endpoint of X (A) * s within 4 units of the corresponding endpoint of <&(B). 
We write A' = x{A)- Let b be one of the endpoints of B. For ease of labelling 
let's assume that A\ is the gene on our list of 75 that has the same type as 
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A. Let Bi be the core of A\ and let b\ be the vertex of B\ corresponding 
to B. Let (3 be the inflation map associated to A\. Let b[ be the endpoint 
of A\ that is within 3 units of B(bi). Let P\ denote the dynamical polygon 
associated to [(Ai,ai)] and let P\ be the dynamical polygon associated to 
[(Ai,&i)]. Note that Pi and P\ are translates of each other. 
Consider a new map f3 defined by the rule 

(3(x) = f3(x + ai — b\) + b[ — a[ 

By Lemmas 13.21 and 13.31 together with the fact that P\ is the relevant trans- 
late of Pi, we see that (3 is a inflation map defined on £(Pi). By the Shadow 
Lemma (3 is 4-pseudo-Lipschitz, and hence so is (3. By translation equiva- 
lence, we have a — b = a\ — b\ and a' — b' = a[ — b[. Therefore 

(3(b) = f3(b + ai- bi) + fe'i - a[ = f3(b + a-b) + b' -a! = b' . (22) 

Note that 

11^)11 = ||6i-*(6!)||< 3 
by the second property of our inflation generator. Therefore 

\\b' - m || =* 0(b) - $(p)|| = ||^(6')|| < ||«^(6) || + ll^(fe') - vtfMW < 

3 + 4 diam(Pi) < 3 + 4(1/4) = 4. 

The starred equality comes from Equation [221 The 1/4 in this last calcula- 
tion comes from the fact that any two points in the same dynamical polygon 
are within 1/4 of each other. A 



3.6 Coherence 

Here we explain how to check that the inflation structure x is coherent. Our 
inflation generator gives rise to a list Pi, P75 of inflation maps, where (3j is 
defined on S(Pj). Here Pj is the dynamical polygon associated to [(Aj, a 3 -)]. 

Say that an extended gene is a polygonal arc of V having length 7. An 
extended gene is just the union of two consecutive genes. Say that an extended 
gene type is an equivalence class, up to translation, of extended genes. It 
turns out that there are 89 extended gene types. Given an extended gene 
X we define xPO = x(Xi) U xP^)) where Xi and X 2 are the two genes 
comprising X. 
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Lemma 3.11 (Coherence) Suppose there is a list Xi,...,X 89 of extended 
genes, representing all the types, such that x{Xj) is a polygonal arc for all 
j = 1, ...,89. Then x is coherent. 

Proof: Let X be an arbitrary extended gene and let X be the translation 
equivalent extended gene on our list of 89. It suffices to show that x(X) and 
XpO are translates of each other. For each object Y we associate to X, we 
let Y be the corresponding object for X. 

Let X\ and X 2 be the two genes comprising X. Let x\ and x 2 be the two 
center points of these genes. Let /3i and /3 2 be the inflation maps associated 
to the gene types of X\ and X 2 . To prove that X is translate equivalent to 
X it suffices to prove 

- P 2 (x 2 ) = AOn) - (3 2 {x 2 ). (23) 

This is what we will do. 

Since \1/ is afline and injective on Z 2 it suffices to show that 

*(/3i(x!)) - V((3 2 (x 2 )) = - *(/3 2 (^))- (24) 

Let pj = ^(xj) and pj = ^(xj). Since X and X are translation equivalent, 
we have X\ — x% = x 2 — x 2 . Since \1/ is afline, we have 

Pi ~ Pi = P2 ~ P2- (25) 

Let 7j : Pj — > T 2 be the special similarity associated to Pj. We work 
in local Euclidean coordinates, so that r yj(x,y) = —<p~ 3 (x,y) + Cj for some 
constant Cj. Then 

7i (Pi) - 7i (Pi) = -0 _3 (Pi - Pi) = -0 _3 (P2 - Pa) = 72(P2) - 72(P2)- 
Hence 

7i(Pi) - 7 2 (p 2 ) = 7i (Pi) - 7 2 (p 2 )- (26) 
But we also know that 

',(/',) = 7i(Pi) = *(&(%))• (27) 

Equations [26] and [27] combine to establish Equation [24] 4 
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4 Computer Aided Verification 



4.1 Overview 

Recall that a gene is a polygonal arc of length 6 contained in To, the com- 
ponent of the arithmetic graph containing (0,0). We are trying to verify 
the existence of an inflation generator. This inflation generator consists of 
a length 75 list of the form {(A,-, A'-)}, where Aj is a gene and A'- is a path 
that nicely shadows Aj in the sense of §3.5. We associate some auxilliary 
objects to our list, namely: 

• Let Pj be the dynamical polygon associated to Aj. 

• Let j3j be the inflation map associated to (Aj,Aj). 

• Let 7j : Pj — > T 2 be the special similarity associated to Pj. Here 
^j o — o Pj, whenever all maps are defined. 

Given these basic objects, here are the things we need to check: 

1. Each endpoint of A'- is within 3 units of the corresponding endpoint 
of &(Bj), the dilation of the core Bj of Aj. We have a list of all the 
vertices involved and we just check this directly. See §8. 

2. For each gene Aj, with center vertex aj, we have correctly computed 
the dynamical polygons Pj associated to [(A,-, a,-)]. 

3. We have r )j{Pj) C Pj, where P'- is the dynamical polygon associated to 
the pointed strand type [(A^,aj)]. Our method will not require us to 
compute Pj explicitly. 

4. There is a complete list X±, X 89 of representatives of extended genes 
such that x(Xj) is a polygonal arc for j = 1, ...,89. Here complete list 
means that every extended gene type is represented. 

5. Our list Ai, A 75 of gene types is exhaustive and our list X\, ...,X 8d 
of extended gene types is exhaustive. There are no other gene types or 
extended gene types. 

In this chapter we will explain the main theoretical points of our verifica- 
tions. In §7 we will write enough pseudo-code so that the interested reader 
can see explicitly how we do all the calculations. 
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4.2 Computing the Polygons 

Let A = Aj denote one of the genes on our list. Let a_ 2 ,...,a 2 denote 
the 5 interior vertices of A. We are interested in computing the dynamical 
polygon P associated to the pointed strand type [(A, do)]. Here we explain 
how P is computed. We don't actually need to know how to compute P for 
our proof— we just need to verify the answer is correct— but the method of 
computation suggests how we verify that the answer is correct. 

We define a dynamical translation of the torus T 2 to be a map of the 
form 

[(x, y)\ - [(:r, y) + (e^" 4 + e^" 3 )]; d, e 2 G {-1, 0, 1}. (28) 

This operation makes sense because T 2 = R 2 / Z 2 is canonically an abelian 
group. We remind the reader that [(x, y)\ is the equivalence class of the point 
(x,y) in T 2 . We say that (61,62) is the type of the dynamical translation. 

Recall that V is partitioned into 26 open polygons (and their boundaries) 
as in Figure 2.5. We associate to A a list Q_ 2 , ■1Q2 °f 5 open polygons of 
the partition V . Here Qj is the polygon having the same local type as the 
vertex aj. Depending on A there are 4 dynamical translations T_ 2 , T_ 1; T 1; T 2 
with the property that x G P if and only if 

• x G Q . 

• 71 (X) G Q x 

• TMTlOr)) G Q 2 . 

• T..i(x) G Q_!. 

• TL 2 (T_i(x)) G Q_ 2 . 

We call these four conditions Property X . Given this information we have 

p = Q n Tf 1 ^ n t 2 \q 2 )) n Tr 1 1 (Q_ 1 n tz 2 1 {q 2 )). (29) 

Billiard King simply computes this intersection. 

If we replace each set Qj by its closure then we arrive at a criterion 
for when x G P. We call this criterion Property X. 

Figure 3.1 (which we repeat as Figure 4.1 below, for convenience) shows 
the 75 dynamical polygons in dark grey. We list the actual coordinates 
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in §8.21 The 24 light polygons comprise the complement. Here is their 
interpretation: A vertex p in the arithmetic graph lies on a closed polygon 
of length 5 or 7 if and only if $f(p) lies in the interior of one of the light 
polygons. We do not use these light polygons in this paper, but we note that 
they "round out" Figure 4.1 by explaining what is left over after we plot our 
75 dynamical polygons. 




Figure 4.1: The partition 



For the first stage of our verification, we use the Billiard King algorithm 
to compute P to high precision. Next, we use the fact that the vertices of P 
are elements of \Z[(f)] to guess the exact expression for the vertices. Let Pq 
denote this guess. Our goal is simply to verify that Pq = P. 

First we check that each vertex of Pq satisfies Property X . This tells us 
that Pq C P. If Pg is a proper subset of P then there is some edge e of Pg 
such that every point x on the interior of e satisfies Property X. We rule 
this out by explicitly choosing one point per each edge of Pg and showing 
that it fails to have Property X. In §7 we explain the calculation in detail. 

When we run the computation for each of the 75 cases it works. Of 
course, the fact that Figure 4.1 is a partition of the torus gives extremely 
strong visual evidence that we have guessed correctly in the first place. 

4.3 Checking the Shadowing Property 

Our method here is quite similar to our method for verifying that Pg = P- 
Let A = Aj be one of our genes. Let (3 — etc. be the objects associated 
to A. Our goal is to check that j(P) C P'. 

Let Q' i, Q' , Q[, ... be the sequence of polygons in the partition V asso- 
ciated to the vertices of A'. We set up the indices so that the Oth vertex of 
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A' is the distinguished point a' that shadows <3>(a). (Here a is the center of 
the gene A.) 

Just as we did for A, we generate the sequence of dynamical translations 
for the strand A'. This is a sequence whose length varies with the choice of 
gene. The length typically varies from 10 to 20. The condition that x G P' 
amounts to checking that x satisfies what we call Property X': 

• x E Q' . 

• T x {x) G Q[ and T^x) G Q'_ x 

• T 2 (T!(x)) G Q' 2 and T_ 2 (T-i(x)) G Q _ 2 

• r 3 (T 2 (T 1 (a;))) G Q' 3 and r_ 3 (T_ 2 (T_ 1 (a;))) G Q _ 3 

• etc. 

We also have the corresponding Property X', the closed version. We 
simply check that each vertex of j(P) satisfies Property X'. This means 
that the closure of j(P) is contained in the closure of P'. Hence "f(P) G P' . 
Billiard King computes all these quantities and displays them visually so that 
the user can see in each case that j(P) C P'. 

4.4 Checking the Coherence 

We will use the forwards direction of To to check the coherence. In order 
to make the construction to follow we first need to know something about a 
certain finite portion of IV We verify by direct inspection that To contains 
a polygonal arc of combinatorial length 2 14 , connecting (0, 0) to a point in 
the positive quadrant. 

In this section we will use the notation and terminology from §3.61 and 
the Coherence Lemma. In particular X±, ...,X$g is a complete list of repre- 
sentatives of the extended gene types. 

For N < 13 let Tq denote the first 2 N segments of T , starting from (0, 0) 
and moving in the direction of the positive quadrant. We call N sufficiently 
large if each extended gene type has a representative on Tq . We show by a 
direct computation that N = 10 is sufficiently large. Once we have our list 
of 89 gene types it is a completely straightforward matter of checking that 
they all occur on Tg A forteriori, the polygonal arc Tj contains each of our 
genes A 1 ,...,A 75 . 
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Our inflation generator gives rise to the inflation structure x- Our goal 
is to prove that x is coherent. Here is the strategy. Given any gene A C Tq , 
we let a be the center point of A. We know that A has the same gene type 
as some gene Aj in our inflation generator. The center point dj of Aj is such 
that $(%) is close to the point (3j(a,j). We record the gene type of the gene 
centered at Pj(aj) and then find a point a' G rg 3 which is near $(a). We 
then verify by direct computation that f3j(a) = a'. 

Now we know that x(A) C IV We check explicitly that the points 
a[, a' 102 4 occur in order on Tq 3 . This means that the two strands x(^-i) 
and x{A-2) overlap whenever A\ and A 2 are consecutive genes on Tq . But 
then x{A\) U x(^) is a polygonal arc. This works for all consecutive genes 
on Tq , including the 89 we need for the Coherence Lemma. Applying the 
Coherence Lemma, we see that x is coherent. See §7 for more details. 



4.5 Checking the Completeness of the Lists 

We compute explicitly that Tq contains 75 gene types and 89 extended gene 
types. We also compute explicitly that Tq 3 has 75 gene types and 89 extended 
gene types. That is, when we go out 8 times as far, we see no new genes 
or gene types. We also recall that Tq is shadowed by a subset of Tq 3 . This 
means that the process of replacing Aj by x(A,), for each of j = 1,...,75 
produces no new gene types and no new extended gene types. Iterating, and 
applying induction we see that our list of gene types and extended gene types 
is complete, in the forwards direction. We then make all the same checks in 
the backwards direction. 

Remark: We could take a different approach to the completeness of our 
list of genes. Figure 4.1 shows a partition of the grey polygons of V into 75 
light grey polygons and 24 dark grey polygons. As we mentioned above, the 
dark grey polygons correspond to points in the arithmetic graph contained 
on closed 5-gons or closed 7-gons. Our list of genes is complete because the 
remaining part of V is completely partitioned by the 75 corresponding dy- 
namical polygons. Any additional gene type would have a dynamical polygon 
that overlaps with one of the ones we already have plotted. We do not insist 
on this approach because it requires an analysis of the 24 dark gray polygons 
in the picture, something we have not attempted. A similar picture would 
reveal the completeness of the list of 89 extended gene types. 
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5 Proof of Theorem 



1.2 



5.1 Existence of the Orbits 

Under our affine map 5" — ► S, discussed in §2.1, the Cantor set C from 
Theorem 11.21 is the Cantor set on the line {y = — 1} with similarity constant 
-3 and endpoints (0, —1) and (20~ 3 , —1). For ease of discussion we identify 
C with a subset of R by dropping the second coordinate. 

Lemma 5.1 The orbit of every point in C* is entirely defined. 

Proof: Lemma 12.11 m the case A = 0~ 3 = 20 — 3, shows that the outer 
billiards map is defined on and preserves (R — 2 Z [<$>]) x Z dd- Hence, it 
suffices to prove that C* n 2Z[(f>] = 0. 

We argue by descent. Suppose x = a + b<p G C*, with a, b G 2Z. We 
must have a ^ 0. We take \a\ > 2 as small as possible. If a < —2 then 

x > = 20~ 3 - x = (-6 - a) + (4 - 6)0 G C* 

has \a'\ < \a\. Hence a > —2. If a = —2 then we must have b > 0. But 
—2 + 20 is already too large to lie in C*. Hence a > 2. 

Case 1: Suppose x lies in the left half of C* . Then 3 x G C*. We compute 

3 x = ( a + 26) + (2a + 36)0. 

By minimality, \a + 26| > a. If a + 26 > a then 6 > and a + 60 > 1. Hence 
x ^LC . If a + 26 < —a then 6 < —a. Hence a + 60 < and x $lC. 

Case 2: Suppose x lies in the right half of C*. The map a(x) = 2 — 3 x 
maps the right half of C* back into C*. We compute 

a(z) = (2 - a - 26) + (2 - 2a - 36)0. 

By minimality |2 — a — 2b\ > a. If 2 — a — 26 > a then a + 6 < 1. Since a 
and 6 are even, this forces a + 6 < 0. But then a + 60 < and x ^ C* . If 
2 - a - 26 < -a then 6 > and a + 60 > 2. Again x ^ C. 

This takes care of all the cases. 4k 
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5.2 A Gap Phenomenon 



Let T_ denote the backwards half of IV We want to see that T_ rises away 
from dH and then back to dH infinitely often (in a precise way). Our main 
idea is to observe that some initial portion of T_ rises and falls, and then 
to propagate this property using the self-similarity. One worry is that the 
approximate nature of the self-similarity of T_ causes the "lowest points of 
approach" to drift away from dH. A gap phenomenon comes to the rescue. 



Equation [7] says that \I/ = ipoT. Recall that P 1; P 2 6 are the open polygons 
comprising the partition V of T 2 . The polygon P 3 is the parallelogram at 
the bottom left of Figure 2.5 and ip([0, 20~ 3 ]) is the long diagonal of P 3 . 

Lemma 5.2 Suppose that x G (0, 16) and ip(x) G P 3 . Then x G (0, 20 -3 ). 

Proof: Looking at Figure 2.5 we see that P 3 is the little parallelogram in the 
bottom left corner. The corner vertex of P 3 is (0, 0) and the opposite vertex 
is (0" 4 ,0- 3 ). So, if ij)(x) G P 3 then 8{x/2) G (O,0~ 3 ). Hence x [20" 3 ,2]. 
It is now an easy exercise to check that the segment ip{[2, 16]) is disjoint from 
P 3 . One way to do this exercise is to plot the ?/>-image of the 2000 maximally 
and evenly spaced points in the interval [2, 16] and observe that none of them 
is within .0001 of P 3 . Since ip is 1-Lipschitz, our plot guarantees that no 
point of [2, 16] lies in P 3 . We omit the details. A 

Given a point p G H, the halfplane containing the arithmetic graph, let 
v (p) denote the length of the vertical line segment connecting p to the line 
through the origin parallel to dH. For instance v(0,m) = m. (The line dH 
lies slightly below the origin.) 

Corollary 5.3 (Gap Phenomenon) Suppose p is a point of the arithmetic 
graph having type 3. If v (p) < 7 then v{p) < 4>~ 3 . 

Proof: Note that T(0, 7) < 15, and the fibers of T are parallel to dH. Hence 
T{p) G (0,15). By definition, 1>(T(p)) G P 3 . Hence T(p) G (O,20~ 3 ). But 
T(0, </>~ 3 ) > 2</r 3 . In other words, if v(p') = 4>~ 3 then T{p') > 2(fi~ 3 . Hence 
v(p) < <p~ 3 . 4 



Let 
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5.3 Infinite Return 



In this section we prove that 0-(p) returns to every neighborhood of p in- 
finitely often. We think of this as a warm-up to the proof of Theorem 11.21 
Let Ao be the gene through the origin. Let A' be the strand that shadows 
it. Figure 5.1 shows the picture. 



50 



Figure 5.1: The gene A Q and its shadow A' . 

The dynamical polygon P associated to A is contained in P 3 , because 
the center of Aq has type 3. Thus Pq is a small tile in the bottom left corner 
of Figure 2.5. Here we show a plot, and describe all the plotted objects 
in order of decreasing size. Each polygon is contained in the previous one. 
Referring to Figure 5.2: 

• The big L is the corner of our fundamental domain for T 2 . 

• The big parallelogram is P3. 

• The big triangle is Pq. 

• The small parallelogram is Pq, the polygon associated to A' Q . 

• The small triangle is 7(-Po)> where 7 is the special similarity fixing 
*(0.0). 




Figure 5.2: The associated dynamical polygons. 
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Let A be some occurrence of A on r_. Let a be the center point of A. 
The first site is (—21,5). Let A' = x(A) be the strand that shadows A. Let 
a' be the center point of A'. Our analysis shows that ^(a) and ^(a') both lie 
in P . The point is that Pq C Pq. Let's assume by induction that v(a) < 0~ 3 . 
Then, by the triangle inequality, and the definition of our inflation structure, 

v(a') < 3 x 0" 3 + 4 = 5. 

Since a' also has type 3, our Corollary above now says that v(a') < -3 . 
We can now start with the copy of A based at a'. This produces a further 
point a", centered on a copy of A , such that v(a") < 0~ 3 . Thus we produce 
distinct points a, a', a",... G T_ all within 0~ 3 of dH. Let a^> be the nth 
point produced by this process and let x^> = T(a^). The points x^ all 
belong to 0_(p). 
By construction 

^(x in+l) ) = 7 (^(x (ri) )); x {n) , x in+1) e (0, 20~ 3 ) C (0, 1/2). 

These two conditions imply that x^ n+1 ^> = ^(x^). Hence x^ = 7 n (x), where 
x = T(— 21, 5). But 7 is a contraction fixing p. Hence x^ converges to p. 



5.4 Cantor Set Structure 



Our argument above came from considering just the gene A . We get the 
Cantor set structure by considering the action of two genes. Notice that A' 
contains two genes having a core of type 3. Figure 5.3 shows the other gene 
A\ with this property. 



V 



Figure 5.3: The gene A\. 

Let A\ be the strand that shadows A\. It turns out that A[ and A' have 
the same type. Moreover, the center point a[ of A[ is the center of another 
copy of Aq. Figure 5.4 shows the corresponding dynamical polygons. This 
time Pi, the little triangle in the bottom left corner, is disjoint from P{. The 
map 7 is the same one as for P , and the parallelogram is again Pq = P[. 
The tiny triangle inside P[ is 7 (-Pi). 
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Figure 5.4: The associated dynamical polygons. 



We set K = P 3 and K = Pq and 70 = 7. There are several geometric 
pieces of information we now record. First, jo(K) = K . Second, let ao and 
di be the two type-3 points of A' . For instance, we might take A' such that 
a = (0,0) and ai = (—5, 1). From Equation [7] we compute that 



V := *(oo) - *(ai) = [(50 



-1 



,50- J )] = (0~ 4 -0 



k-7 a-3 



(30) 



Let t{x) — x — Vq. Geometrically, the translated parallelogram K\ = t(Kq) 
has (0, 0) as a vertex and fits exactly into the lower left corner of K. Figure 
5.5 shows a schematic picture. Let 71 = ro^g. Then jj(K) = Kj. We define 
^00 = lo(K ) and K 01 = 70(^1) 
of the first few of these sets. 



etc. Figure 5.5 shows a schematic picture 



K 



Kll 



K10 



KO 



KOI 



K00 



Kl 



Figure 5.5: A schematic picture 

The reader might expect Kqq to be on the top and K\\ on the bottom, 
rather than in the middle. This twisting of the labelling comes from the fact 
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that the multiplier of 7 is — -3 , a negative number. The fixed point of 7 is 
the nested intersection K n K 00 fl i^ooo--- The limit set C of the semigroup 
(7o,7i) is a Cantor set whose endpoints are the far opposite vertices of K . 

Now we can explain the Cantor set structure. Suppose that X is a type 
3 point T_ such that X is the center of a gene A that is a copy of either Aq 
or A 1 . Suppose also that v(X) < (f)~ 3 . Let X and Xi be the two points 
of A' having type 3. Our notation is such that X Q is the center of a copy 
of the gene A and X\ is the center of the copy of the gene A\. The point 
X shadows the dilated point (p 3 X, and the point X\ is several units to the 
right. Our arguments above show that 

• v(X ) and v(X\) are both less than 0~ 3 . 

• #(Xt) C Ki and *(X ) C if - 

• r(*(X )) = 

. *(X i )= 7i (*(X)). 

Now we consider a tree-like induction process. We start with the point 
X = (0,0) G r_. Considering the above inflation process, we produce the 
points 

X =X; X 1 = (-5,l). 
Inflating again, we produce the points 

-^00 — Xq\ Xoi = X\\ Xxo = (—21,5); Xu = (—26,5). 

Inflating again, we produce the 8 points, corresponding to the binary strings 
of length 3, and so on. In general, we produce one point Xp G T_ for each 
finite binary string (3. (These points are not all distinct, e.g. X 00 = X .) By 
construction 

*{Xfi) C K p . 

Therefore, the closure of the vertex set of ^(r_) contains the Cantor set C . 

Let Iq = [0, 20~ 3 ]. We have already remarked that i/j(Io) is long diagonal 
of P3, the segment that connects the two endpoints of C . Recall that T 
maps the vertices of T_ into the backwards orbit 0-(p). Also f = f o T. 
Therefore C = ip{C), where C C Jo is an affine image of C. Indeed, C is 
the Cantor set from Theorem 11.21 (We are identifying [0, 00) x {—1} with 
[0, 00), as discussed at the beginning of the chapter.) We have realized C as 
a set contained in the closure of Hence 0_(p) returns densely to C. 
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5.5 2-adic Structure 



Since ip(Io) is the long diagonal of P3, the two conditions x E Z 2 n H and 
T(x) < 2<j)~ 3 imply that x has type 3. We say that such points are basepoints 

of r. 

We see inductively that our shadowing construction accounts for all the 
basepoints of T_. That is, every such point arises on a strand of T_ that 
shadows an inflated gene. Therefore, we can index the basepoints by binary 
strings (3. Each binary string {3 represents the integer n{(3) as usual. For 
instance 

n(11001) = (1,1,0,0,1) • (16,8,4,2,1) = 25. 

Examining our construction we see that the two basepoints Xp and Xp/ are 
equal if and only if n(/3) = n(/3'). Moreover, as we travel along T_ away from 
(0, 0) we encounter the points Xp in order of the integers they represent! 

Let 62 '■ Z2 — > C be the homeomorphism from Theorem II .21 By construc- 
tion, we have ^ 1 (0) = p. Let T : Z 2 — > [0, 00) be the map from Equation [3j 
Let X and X and Xi be the points referred to in our shadowing construction 
above. We have already mentioned that 

*(x )=7iP0; j = Q > 1 - 

Given the defining property of 82, we have 

e^tnXj)) = 29^(T(X))+j; j = 0,1. 
But then we have 

e 2 {T{Xp))=n{(3). 

But, by construction T(Xp) is the nth point of 0-(p) fl/o, which is the same 
as the nth point of 0_ (p) fl C. Hence 6* 2 maps the nth point of 0_ (p) n C 
to n. This is a special case of the first statement of Theorem 11.21 

Now we will deal with the question of return times and excursion distances 
for the points of 0_(n). Given a basepoint X, let T + (X) be the forwards 
portion of To which starts at X. Let X and Xq and X\ be as in the previous 
section, so that Xq is the basepoint that shadows <p 3 X. 

By the inflation property, r + (X ) closely follows (f) 3 T + (X). Also, one can 
see by looking at a single example that the initial portion of r + (X ) rises up 
at least 10 units from OH before coming back towards OH. Moreover, the 
next basepoint (after Xq) encountered by r + (X ) is the one that shadows 
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the dilation of the next basepoint encountered by F + (X). These properties 
show that r + (A ) rises up about 3 times as high as T + (X), and takes about 
<f) 3 times as long to return to the next basepoint as does T + (X). The next 
basepoint encountered by r + (A ) is the one that shadows the dilation of 
the next basepoint encountered by T + (X). On the other hand, r + (Ai) just 
travels a few units to the right before returning to the basepoint T + (X). 

Given any binary string (3, let z/(/?) denote the number of Os on the right 
of j3. For instance z/(11000) = 3. Equivalently, z/(/5) equals the highest power 
of 2 dividing n{[3). Given X = Xp let v = vp. Applied inductively, our 
arguments show that T + (X) rises up roughly (p 3u units before returning to 
the next basepoint after roughly (j) 3u units of time. Here X = Xp. Hence the 
excursion distances and return times for the forward orbit of the nth point 
x n of 0-(p) fl C are proper functions of the 2-adic distance from 6 , ^" 1 (a;„) to 
0. 

By simply reversing the direction of T + (X) we see that the excursion 
distances and return times for the forwards orbit of x n are proper functions 
of the 2-adic distance from 6 l 2 (x n ) to —1. The point here is that 62 maps n 
and n — 1 respectively to the endpoint of T + (j3) and the first basepoint it 
encounters when travelling to the right. 

There is one more observation we want to make. If X is a basepoint and 
v{X) is very large, then a very long initial portion of T + (X) looks exactly 
like T + (0, 0), the forward portion of T that starts at (0, 0). This is because 
both strands are "produced" by many iterates of the same inflation process. 
An equivalent way to see this is that any finite portion of To is determined by 
some small dynamical polygon in T 2 containing \l/(0,0). If X is some point 
with v(X) large, then ^(A) belongs to this same dynamical polygon. 

5.6 Extension to the Cantor Set 

Now we put everything together and finish the proof of Theorem 11.21 

Lemma 5.4 (Rising) For each positive constant K and each lattice point 
X G there exists a constant K' = K'(K,X) with the following property: 
If we take K' steps in either direction along Y_, starting at X, then we rise 
at least K units above OH . 

Proof: Consider the forwards direction. The backwards direction is similar. 
It this lemma is false then we can encounter a long string of basepoints {Aj} 
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such that v{Xi) is always small. But the map x — > x — 1, when iterated, 
brings any point very close 2-adically to within a uniform number of steps. 
Hence, we don't have to walk very long from X before we hit a basepoint 
with high z/-value, and then we rise up steadily away from dH for a long time. 
This comes from the fact that a long initial portion of T + (X) is a translate 
of a long initial portion of T + (0, 0), as we remarked above. 4k 

Lemma 5.5 For any y G C* and any K > there are positive integers 
n + = n + (y,K) and n_ = n^(y,K) such that the (n + )th point of + (k) and 
the (ri-)th point of 0-(y) are at both least K units from the origin. 

Proof: We can find a sequence {x n } G 0-(p) converging to y. For n suffi- 
ciently large, the first K' iterates of x n and y will have the same combinatorial 
structure. Here K' — > oo as n —>■ oo. But then the first K' iterates of y re- 
main (say) within 1-unit of the corresponding K' iterates of x n . Letting K' 
be as in the Rising Lemma, we see then that some point of x n rises up K 
units and therefore some iterate of y rises up at least K — 1 units. Since K 
is arbitrary, the orbit of y (in either direction) is unbounded. A 

We continue with the notation from the lemma. Let i n _ denote the 
forward return to C* of the point x n . As long as y ^ p, there is some uniform 
K such that x n and separated by at most K units. Here K depends 

on y but not on n. This we see that y_ and y are separated by at most K 
units. For n large enough the orbit of x n has the same combinatorial structure 
as the orbit of y for the first K iterates. Hence y^—y = x n _ —x n . Therefore, 
by continuity, 0% = 9^ (y) — 1. The same continuity argument takes 

care of the statement about the return times and excursion distances for the 
forward orbit. This proves the first statement of Theorem II .2\ and the second 
statement has essentially the same proof. 
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6 Proof of the Arithmetic Graph Lemma 



6.1 The Computational Evidence 

Before we launch into the proof of the Arithmetic Graph Lemma we discuss 
how we discovered that it was true. Our proof of the Arithmetic Graph 
Lemma is logically independent from the way we discovered it, but the com- 
putations we made during the discovery process will serve us below at a 
certain step in our proof, namely the step where we show that Lemmas 16.61 
and 16.71 below are equivalent to the Arithmetic Graph Lemma. 

Recall that r C H, a certain halfplane in R 2 . We thought of T as a 
union of paths of physical particles, and we got the idea to partition H into 
parallel bands, each of which has width 0. We thought of these bands as semi- 
permeable membranes, letting various particles pass through and deflecting 
others. Experimentally, we noticed that the permeability of the nth band 
depended somehow on the decimal part of n/<p. In particular, the bands 
indexed by Fibonacci numbers are the least permeable and behave the most 
like mirrors. Thus, the quantity 5(n/<p), an R/Z valued coordinate, seemed 
relevant to the "physics" of a particle in the nth band. Here 5(x) is the 
decimal part of x. 

The other quantity that seemed to influence the "physics" of a particle 
was the relative position of the particle within the band that contained it. 
Like our first coordinate, this is another coordinate that takes values in R/ Z. 
Thus, the local properties of T seemed to be determined by two coordinates 
in R/Z, which we considered as a point in T 2 . When we wrote down the 
formulas for these coordinates and tuned them, we arrived at our map \l/ 
from Equation [71 

Once we had \l/ we performed an experiment. We took some huge sample 
of T, and fixed some local type. Sampling lattice points p, we plotted ^(p) if 
and only if p had that local type. We noticed that the resulting plot was very 
dense in either a single convex polygon, or else the union of several convex 
polygons, depending on the type. (Later we indexed the types so that each 
type corresponds to a single polygon.) It was easy to guess the vertices of 
the polygons from the plots, and this is how we determined the partition V. 
Indeed, Billiard King draws T in two ways, using the partition and using the 
dynamics, and the picture is the same. Now we turn to the proof. 

5 The directness of this discussion is misleading. We arrived at the ideas here only after 
having exhausted every half-baked and useless scheme we could think up. 
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6.2 Factoring the Return Map 

Say that a strip is a region A C R 2 bounded by 2 parallel lines, 8qA and 
d\A. Let V be a vector such that 8qA + V = d\A. Given the pair (A, V) we 
(generically) define a map E : R 2 — ► A by the formula -E(x) = x — where 
n is the unique integer such that E(x) G A. This map is well defined unless 
x lies in a discrete infinite family of parallel lines. 

There is a unique affine functional f(x,y) = a%x + a^y + 03 such that 
h{V) = 1, and f(x,y) G (0, 1) iff (x,y) G A. Here f L {x,y) = a x x + a 2 y is 
the linear part of /. Given / we have the following explicit formula: 

E(p) = p - floor (/(p)) V. (31) 

Equation ED is defined unless f(p) is an integer. We say that a = (ai,a 2 ,a 3 ) 
is the triple associated to (A, V). 




Figure 6.1: The pinwheel 

Each strip Aj in Figure 6.1 is obtained by extending an edge of the side 
of S, and then rotating this extended side through one of the vertices that 
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does not contain the edge. Each vector Vj is twice the difference between 
a pair of vertices on the kite. Let Ei, E 2 , E$, E± be the corresponding strip 
maps. Here is the data associated to these maps: 



Let E j+A = Ej for j = 1, 2, 3, 4. 

Lemma 6.1 (Pinwheel) Let x E C(±). Let a^- = Ej(xj-i). Then x 8 and 
Yr(xo) on the same vertical line. 

Proof: Let Vj+4 = —Vj for j = 1,2,3,4. For any x on which T is defined, 
we have T(x) — x = V x , where V x is twice one of the vectors pointing from 
one vertex of S to another. The vectors V±, Vg have this form. 

The lines comprising our strips divide R 2 into a finite number of bounded 
regions and a finite number of unbounded regions. Let K denote the closure 
of the union of bounded regions. One can check, with a little experimentation, 
that the finitely many iterates p, T(p), Yr(p) avoid the set K provided that 
p = (xq, ±1) and x > 4. We will consider such points first. 

The complement of K is divided into 8 sectors Si,..., Eg. Each sector 
Ej is bounded by one line of Aj_i and one line of Aj, and contains a non- 
compact subset of Aj-i. The left hand side of Figure 6.1 indicates half of 
these sectors. It is easy to check that q G Ej implies V q = Vj. This result 
immediately implies our lemma for points (x , ±1) where x > 4 and all maps 
are defined. 

When xq < 4 we can just check the few cases by hand. Indeed, to be sure 
we sampled 1 million evenly spaced points in the two intervals (0, 10) x {±1}. 



The Pinwheel Lemma lets us factor the return map. The composition 
E S ...E 1 maps x E C(±) to x s E R 2 , and the ^-coordinate of x$ is ±1 + Ak. 
In each case, ^r(xq) E C(±). Hence 



a 1 = (-1/4, +1/4, +3/4) 
a 2 = (-0/4, +1/2 -0/4, +1/2 
a 3 = (-0/4, -1/2 -0/4, +1/2 
a 4 = (-1/4, -1/4, +3/4) 



0/4) 
0/4) 



V 

V 2 

V 3 

v, 



(0,4) 
(-2, +2) 
(+4-40,0). 
(-2,-2). 



T ii = (o(£ 4 o...o J B 1 ) 2 ; 



C(x,±l + 4A;) = (x,±l). 



(32) 
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6.3 Four Dimensional Compactificaton 

We say that a strip map E, with associated triple (01,02,03) and vector 
V = (v\,v 2 ), is special if 

• ai, a 2 , a 3 G and a\ ± a 2 G \Z[<f>\. 

• Vi, V 2 G 2Z[0] and ^ ± w 2 e 4Z [</»]. 

By inspection, the strip maps from Equation [32] are all special. 

Let Tg = (R/8Z) 4 . Consider the following embedding of R 2 into Tg 4 : 

^(^ 2/) = [(a; + y, ar - y, (2 + y)/<j>,(x- y)/4>)]%- (33) 

Here [i>]g denotes the image of v G -R 4 in Tg. We only care about the next 
result for our 4 maps, but the argument is easy to make in general. We give 
explicit formulas for the extensions of our 4 maps in §7.81 

Lemma 6.2 ( Extension) Let E be a special strip map. There is a finite 
union Y of flat 3 -dimensional tori in Tg together with a map E : Tg — Y — > Tg 
such that ip o E = E o ip whenever both maps are defined. E is locally affine 
on each component ofTf — Y and the linear part of E is defined over Z[(f>] 
and independent of component. 



Proof: We fix p\ G R 2 on which E is defined, and consider points p 2 such 
that ip{j>\) and ip{j>2) are close together in Tg . Setting p = p 2 —pi, this means 
that the coordinates of ip(jp) are all near 0. We make this assumption about 
the coordinates of ip(p) and treat these coordinates as small real numbers. 
Let p = (x, y). Given Tj G Z for j = 1, 2, 3, 4 we write 

<— )-ra-(^)^(^)^(^) <*> 

The coefficients of Tj are nearly integers. Hence S((ti, ...,t 4 )) is near 0, and 
also is a |Z-linear combination of the coordinates of ip(p). We call this the 
reduction principle. Also, <f>{ri, r 2 , r 3 , r 4 ) = (r{, r 2 , r^, 74) because = 
Call this the multiplication principle. 

Let / be the functional associated to E. Since Z[<p] = Z[(f)~ l ] we have 

11 11 
ai = -s 1 + —ti ) a 2 = -s 2 + —t 2 ; s±, t±, s 2 , t 2 E Z. 
A Am A Ad) 
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The special property of our map implies that u± = \(u\ ± u 2 ) G Z for 
u = s,t. We have 

f( P 2) - f(pi) = (si I + Sa |) + (ti^ + t 2 ^) = 2( S+ , a _, f + ,f_). (35) 

From equation [35] and the reduction principle, f(pi) — /(P2) is extremely near 
an integer, and much closer to an integer than /(pi). Hence /(P2) is n °t an 
integer and is defined on p 2 . 

Let iV 3 - be the floor of f(pj)- The properties of /(pi) — /(P2) and /(p?) 
just mentioned imply that N 2 — Ni is the integer closest to f(p 2 ) — /(pi)- 
Equation [351 gives us N 2 — Ni = 2K with 

K = (3-5((3)eZ ] p=(s + ,s.,t + ,t_). (36) 

Let V = (vi,v 2 ) be the vector associated to E. Equation [311 yields 

E{p 2 ) = E( Pl )+p-2KV. (37) 

To prove the existence of E it suffices to show that each coordinate of ip(2KV) 
is a Z[(p] linear combination of the coordinates of ip(p). Let 2KV = (wi, w 2 ). 
By hypotheses Vi±v 2 G 4Z[0]. Hence (p a (wi±w 2 ) G 8Z[0] for any exponent 
a G Z. (We only care about a = 0, — 1.) In light of this observation, it 
suffices to show that S(K(p) is near and also a linear combination of 

the coordinates of ip{j>). 
Equation [36] gives 

8{K<t>) = 8(W - MP)) = Sm - 5(08(13)) =• 5m - <j>5((3). (38) 

The starred inequality comes from the fact that 5(/3) is near 0. From the 
reduction principle and the real multiplication principle we see that both 
5(j3) and 5(<pj3) are near and ^Z linear combinations of the coordinates 
i/j(p). This proves what we want about the coordinates if %jj(2KV). 

It only remains to find the domain of definition for E on T 8 4 . The map 
E is defined on the complement of an infinite union of evenly spaced lines in 
R 2 . These lines consist of points (x, y) satisfying the integral linear equations 
(s 3 , S4, t 3 , t 4 ) G Z. Applying the map ip we see that the image of these points 
is contained in a finite union of flat embedded 3-tori. Our construction shows 
that E is defined and locally affine on each component of Tg — Y. 4t 
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6.4 Geometry of the Return Map 

We begin with a consequence of the Pinwheel Lemma. Let iti : R 2 — ► R be 
projection onto the first coordinate. 

Lemma 6.3 The map tci(Tr(x, 1)) — x only takes on finitely many values. 

Proof: To see this note that the "octagonal spiral" S in Figure 6.1 connect- 
ing the points xo, Xg, C(^s) is within a bounded distance of being symmetric 
about the origin. In other words, the image of E under reflection through the 
origin lies within a uniformly bounded tubular neighborhood of S. In terms of 
the displacement vectors X\— x , x 2 —xi, we see that the opposite sides of £ 
cancel, up to a uniformly bounded error. Hence there are uniformly bounded 
integers ctj, depending on (x, 1), such that Tr(x, 1) — (x, 1) = J2 a jVj- 4 

We now return to the theme taken up in the previous section. The analog 
of the Extension Lemma is not quite true for the map ( in Equation [32] 
but nonetheless a similar result holds. The domain for ( is the union L 
of horizontal lines whose y-coordinates are odd integers. Consider the map 
if) : C(±) -> T 2 given by 

ifj(x,±l) = [(x/(2 ( j ) ),x/2)} (39) 
Lemma 6.4 There is a locally affine map ( : Tg — > T 2 such that (oijj = ipo(. 

Proof: We define ( : R 4 — > T 2 by the formula 

~ \x 3 + x 4 Xi + x 2 

((X 1 ,X2,X 3 ,X 4 ) = , 

The image in T 2 is unchanged if we add multiples of 8 to the X; coordinates. 
Hence ( factors through a locally affine map (which we give the same name) 
C : Tg — » T 2 . Letting m be an odd integer, we compute 

~ ~ ~( x + m x - m\ 

Q o ip(x, m) = Q\x + m, x — m, — — — , — — — j = 

[(x/2(j)),x/2)} = ^{x, ±l)=i)o ((x, m). 
The choice of sign for ±1 depends on the congruence of m mod 4. 4jk 
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Let 

T 2 (±) = closure(^(C(±))) C T 8 4 . (40) 

Notice that the first and third (or second and fourth) coordinates of ip com- 
pletely determine the whole image of C(+) in Tg. Thus T 2 (+) is a flat 
2-dimensional torus. The same goes for T 2 (— ). 

Let Ej be the extension to Tg of Ej. We give the explicit formula in §7.81 
Referring to Lemma 16.41 we define 

f R : f 2 (±) - T 2 ; f R = ( o (E 4 o ... o E^ 2 . (41) 

Combining the Extension Lemma with Lemma 16.41 we have 

T R o^ = ^oT K (42) 

Lemma 6.5 (Constancy Lemma) Suppose that X C T 2 (±) is a path 
connected open set on which the map T R is entirely defined. Suppose that 
Pi,P2 £ C(±) are points such that ip(pj) £ X for j = 1,2. Then 

tti(Tr(pi) -pi) = HiiJ R {p 2 ) -p 2 )- 

Proof: Assume without loss of generality that the path is on T 2 (+). Since 
T R is entirely defined on X it is also continous on X. The map 

v(a) = f R (a)-((a)eT 2 

is a continuous function of a £ X. 

A dense subset of X has the form ^(?/, 1), where y E R. For such points 
we have 

v$(y, 1)) = x T fl o 1) - C o $(y, l)= 2 *Po T R (y, 1) - ^ o C(y, 1) = 3 

o T R (y, 1) - 1) = 4 ^ o n(T R (y, 1) - y, 1) (43) 

The first equality is by definition; the second one is Equation H2] and Lemma 
I6.4t the third one is the fact that ( is the identity on C(±); the fourth one 
comes from linearity, and from the fact that ip only depends on the first 
coordinate. By Lemma I6.3[ the last expression in Equation H3] takes on only 
finitely many values. Therefore, v (y, 1) takes on only finitely many values. 

Now we know that the continuous map v only takes finitely many values 
on a dense subset of the path connected subset X. Therefore v is constant 
on X. Since v is constant on X, the calculation in Equation H3] gives the 
same answer for all y £ C(+) such that ip(y) £ X. The conclusion of the 
lemma follows immediately from setting y — X\, x-i. 4k 
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6.5 The Proof Modulo the Question of Definedness 

Referring to Equations [7] and [391 we have — ip o T, assuming that we 
identify the range of T with either C(+) or C(— ) by padding the second 
coordinate with a ±1. Given x + G let x_ be the point in C(— ) with 

the same first coordinate. Let 

v{x±) = 7n(Tfl(x±) - x±). (44) 

Lemma 6.6 The point ijj(x + ) = ip(xJ) lies in an open polygon of the parti- 
tionV ofT , and this polygon determines the unordered pair {v(x + ),v(x )} . 

Computing just a single point per polygon— we computed millions— we 
check that the quantities given in Lemma [6.61 when hit with the map T , 
match the types shown in Figure 2.5. Therefore Lemma [6.61 (and this small 
computation) implies the Arithmetic Graph Lemma. 

The polygons of V do not determine the ordered pair of numbers in 
Lemma 16.61 For this we need to pass to a finite cover. The map C, : Tg — > T 2 
from Lemma [6.41 gives a finite covering map from T 2 (+) to T 2 . We lift the 
partition V to a finite partition P of T 2 (+). We do the same thing for T 2 (— ). 

Lemma 6.7 The point ip{x + ) lies in an open polygon ofV and this polygon 
determines v(x + ) . Likewise for v (x_) . 

A small amount of computation shows that Lemma 16.71 implies Lemma 
16.61 Hence Lemma 16.71 implies the Arithmetic Graph Lemma. 

Lemma 6.8 Suppose that T^j is defined on each open polygon of the partition 
V o/T 2 (±). Then Lemma \6. 7| is true. 



Proof: It follows immediately from Lemma 16.81 and from the Constancy 
Lemma that if ip(x + ) lies in an open polygon of V, then v(x + ) is determined 
by this polygon. Likewise for x_. We just have to rule out the possibility 
that ip(x + ) lies in the boundary of one of the polygons. The point x_ has 
the same treatment. 

Suppose that ip(x + ) lies in the boundary of some tile of V . We know 
T R is defined at x + . But then T R would be defined and continuous in a 
neighborhood of ip(x + ) on T 2 (+). But then two adjacent tiles of P on T 2 (+) 
would determine the same v values. But then two adjacent tiles of P would 
determine the same unordered pair of v values. We check computationally 
that this does not happen. Contradiction. 4jk 
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6.6 Verifying the Definedness 

At this point we have reduced the Arithmetic Graph Lemma to verifying the 
hupotheses of Lemma 16.81 We will explain how we make the verification. 

We know that Ej+i = Ej for j = 1, 2, 3, 4, but for notational purposes, it 
is convenient to just write E ly ...,E 8 . We let F k = E k ...Ei. We interpret F 
as the identity map. Also, T R = ( o F 8 . To show that T R is defined at some 
point x it suffices to show that the maps F\, F$ are all defined on x. 

We will see from the explicit description of our map E k , given in §7.81 that 
there is a union Y k of two 3-tori, such that Tg — Y k consists of two connected 
components C k (0) and C k (l) on which the map E k is entirely defined and 
locally affine. 

If is defined on a point j 6 T 8 4 then we can define a length 8 binary 
sequence e\, ...,€$, by the property that F k -i(x) G C k (e k ) for k — 1, ...,8. In 
short, we can associate a canonical binary sequence of length 8 to any point 
x on which T R is defined. We call this sequence the itinerary of x. 

Here we define a slightly more general notion of an itinerary. Given an 
itinerary ei,...,e§ we can define E k on Y k by demanding that E k extends 
continuously to the closure of C k (e k ). Then E k is completely defined on Tg 4 . 
It is continuous on the closure of Cfc(efc) and on the interior of (7^(1 — e^), 
but not globally continuous. With this definition, we say that x G Tg has 
extended itinerary ei, e 8 if the maps F k (when extended) are all defined on 
x and F k _i(x) lies in the closure of C k (e k ) for k = 1, 8. Put another way, 
x has extended itinerary e provided there exist points arbitrarily close to x 
that have itinerery e. 

We define the stretch of a convex polygon P to be the maximum distance 
in Tg between consecutive vertices of P. We denote this by c(P). 

Lemma 6.9 (Definedness Criterion) Suppose that 

• Tr is defined on some point of P, and this point has itinerary e. 

• All the vertices of P have extended itinerary e. 

• The stretch of F k (P) is less than 2\/2 for each k = 0, 7. 

Then T R is defined on all points of P, and all these points have itinerary e. 

Proof: Looking at explicit formulas for our maps given in §7.8.11 we see 
that the Euclidean distance between separare components of Y k is at least 
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2y/2. Hence, if Vi,V2 are two points in the closure of Ck(ek) which are less 
than 2y/2 apart, then the line segment W[V2 lies in the closure Ck{^k)- 

Suppose we have shown by induction that Fk-i is defined on P. Let 
P' = Fk-i(P). The hypotheses of this lemma say that every edge of P' has 
length less than 2 a/2, and the endpoints of such an edge are in Cfc(efc). Hence, 
all the edges of P' lie in Cfc(efc). Since dCk(ek) is three dimensional, and (by 
induction) P' is a planar polygon, we must have P' C Ck{£k)- 

If some point of the open P 1 actually lies in dCk{tk) then all of P' must 
lie in dCk{^k)- The point here is that the tangent plane at this bad point 
must be contained in the tangent space to Y^, because there is no crossing 
allowed. We also know that some point of P' lies in Cfc(efc), so the above bad 
situation cannot occur. Hence all points of P' lie in C^ej.). This completes 
the induction step. 4 

We now mention a trick that makes the Definedness Criterion more useful. 
Given a line segment s £ T 8 4 , with endpoints p\ and P2, we let s' denote the 
partition of s into the two segments \pi,q] and [q,Po\ where 

q = Pi<p~ 2 +p 2 <f>~ 1 - (45) 

Note that q G s because _1 + <p~ 2 = 1. We might have chosen q to be the 
midpoint of s, but our choice interacts better with Z[<f>]. Given a polygon 
P, we let P' denote the polygon, with twice as many vertices, obtained by 
subdividing each edge of P. In general, let P^ n > = (P^ 1 ^)'. Then ?W 
has 2 n times as many vertices as P and the stretch of P^ is _n times the 
stretch of P. Thus, by taking a sufficiently large integer n, we can guarantee 
the last condition of the Discreteness Criterion without even computing the 
stretch. 

We take each polygon P of our partition V and perform the following 
calculation. We take the 10th subdivision p( 10 ) and check the Definedness 
Criterion. (This is overkill.) This verifies the hypotheses of Lemma [6.81 and 
thereby completes the proof of the Arithmetic Graph Lemma. 

We give details of our calculation in §7.81 
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7 The Code in Detail 



7.1 Arithmetic Operations 

We do all our computations with complex numbers of the form 

Xo + X 1 (f) X 2 +X 3 (f) 

—r~ + 1 —2— (46) 

Here I = \J — 1. We call such a number an IntegerComplex and represent it 
as a sequence 

{x ,xi,x 2 ,x 3 }. (47) 

Given a and b, both IntegerComplex objects, we perform the following op- 
erations: 

• a + b and a — b are computed by adding or subtracting components. 

• 2ab, another IntegerComplex, is computed by expanding out all the 
terms and grouping them. 

• The conjugate a is obtained by negating the third and fourth coordi- 
nates of a. 

• We have a routine interpolate^, o), which computes a0 _1 + b(fi~ 2 . 

Now we explain how we test the sign of the number a = a + ai0. Since 
(f> is irrational, we need a trick. The idea is to consider 

/so = 12586269025; f 51 = 20365011074 / 52 = 32951280099. (48) 

Here f n is the nth Fibonacci number. The quantity a + ai<p is positive 
(respectively negative) provided that both sums 

si(a) = ao/50 + ai/51; s 2 (a) = ai/51 + a 2 / 5 2 (49) 

are positive (respectively negative). This works because the sign of the differ- 
ence — f n +i/f n alternates with n. We implement our routine, called sign, 
using the Biglnteger class in Java, which does integer arithmetic correctly 
on huge integers. The routine sign certainly could fail for some inputs, but 
it never fails for the inputs we give it. The point is that we always give it 
inputs involving pretty small integers. 
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Our routine dec takes the decimal part of a Integer-Complex and returns 
a Integer-Complex in [—1/2, 1/2]. In defining dec we take the decimal part 
of a IntegerComplex using floating point arithmetic to compute the nearest 
integer. The computation does not use exact integer arithmetic, but then 
we use the sign routine to check rigorously that our guess always lies in 
[—1/2,1/2], thereby guaranteeing that get the same answer as if we had 
used purely integer arithmetic calculations. 

With these preliminaries in place, we describe the IntegerComplex ver- 
sion of our map ^ given in Equation [7J 

psi(a,fo): 

let x = 4 — 12a + 46 and y = 8a — 2. (the IntegerComplex version of Eq[3]) 
Let ci = {(—x + y)/2, x/2}. (division by 20; works because x, y are even) 
Let C2 = {x/2,y/2}. (division by 2; works because x, y are even) 
Let dj =dec(cj) for j = 1, 2. 
return(di + Id 2 ). 



7.2 Classification into Types 

An IntegerPolygon is a finite list of IntegerComplexes, namely the ver- 
tices. Given an IntegerComplex z and an IntegerPolygon P, let Ti(z,P) 
denote the triangle determined by the ordered triple of IntegerComplexes z, 
P(i) and P(i + 1). The indices are taken cyclically. Using our sign routine, 
we have a straightforward routine signArea that computes the orientation 
of a triangle of IntegerComplexes and vanishes of the points are collinear. 
The following routine returns a 1 if z is contained in the interior of P and a 
otherwise. A straightforward variant, IsContainedClosed checks if z £ P. 

isContainedOpen(^,P) : 

loop for the number of vertices of P: 

check that Ti(z,P) and T i+ i(z,P) have nonzero signArea 

if(false) return(O) 

check that Ti(z,P) and Tj + i(z,P) have the same signArea 

if(false) return(O) 
endloop 
return(l) 
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Recall that the arithmetic graph T is a certain subset of Z 2 D H, where H 
is a certain half-plane. Each point (x, y) G Z 2 has a type, which we compute 
by determining which open polygon of V contains ty(x,y) G T 2 . Here we 
explain how we compute this in practice. Lifting V to R 2 we get a tiling V 
of R 2 . We treat psi(x,?/) as a point in i? 2 and check which tile of V it lands 
in. We have arranged that the image of psi is fairly close to the origin, and 
so we only have to check a smallish portion of the tiling. 

We accomplish our goal using two routines. The first of our routines 
checks whether or not an IntegerComplex z near the origin in R 2 is con- 
tained in a union of integral translates of a given IntegerPolygon P. The 
program returns a 1 if the answer is yes. 

isLatticeContained0pen(,2,P): 

loop over i from —3 to 3 and over j from —3 to 3 

if(isContained0pen(;2 + (i,j),P)= 1 then return(l) 
endloop 

We always take P as one of the polygons from our partition V. These 
polygons are listed in §7.51 (See also Figure 2.5.) Our routine classify com- 
bines our routine psi with the routine isLatticeContainedOpen to classify 
each point of Z 2 into the local types. 

We use the routine isContainedClosed in place of isContainedOpen in 
case we want to check that a given point is contained in a closed polygon. 

7.3 Types and Dynamical Translations 

There are 23 nontrivial local types of vertex in the arithmetic graph, as shown 
in Figure 2.4. We list these types here (rather than in the appendix.) The 
array 

T 

oi h (50) 
a 2 b 2 

indicates that one of the edges emanating from a vertex of type Tj. is (ax, o, 2 ) 
and the other one is (&x, b 2 ). In other words, the two columns of the matrix 
encode the type. The ordering of the two columns is arbitrary. We list 
the types T 1; T 3 , ...,T 23 . Similar to what happens above, T 2 j is obtained by 
negating all the entries of T 2j _x, for j = 1, 11. Here is the list: 
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T-i 




T 3 




T 5 


T 7 




1 





-1 


1 


-1 


-1 





1 


1 


1 


1 


1 -1 


1 


-1 


To 




T u 






T 15 







-1 


1 


-1 


1 -1 





-1 


1 


-1 


1 





1 


-1 





Tvt 




Tig 




T 21 


T23 







-1 





-1 











-1 





1 





1 -1 


1 


-1 



For each of these types, there are two associated dynamical translations 
of T 2 . Referring to the matrix in Equation the two maps are: 

0, y) - [( ai 0- 4 + a 2 ct>-\ ai 0- 3 )]; (x, y) - [(b^ + 6 2 0~\ &10" 3 )]- 

(51) 

These maps relate the points ^>(v) and \I/(V) where w is the vertex and v' 
is one of the vertices connected to v by the arithmetic graph. The following 
routine starts with a pair (x, y) obtained from one of the columns of the 
above matrices and returns the IntegerComplex z which effects the corre- 
sponding dynamical translation. That is, on R 2 the dynamical translation is 
given by w — > w + z. The next routine gets the IntegerComplex by which 
we translate. Here x, y G {—1, 0, 1}. 

getMap(x,y): 

return the integer complex with coordinates (lOx — 2y, —6x + 2y, —6x, 4x) 
7.4 The Sequence Generator 

If v 1 and f 2 are two consecutive vertices of the arithmetic graph, then one of 
the two maps associated to v\ coincides with one of the two maps associated 
to t>2- For instance, if a vertex of type 1 is connected to a vertex of type 21 
we could write 

n t 21 1 21 

10 or 10 

11-1 11 

The second notation system is a simplification of the first. We drop off the 
last column because it is not of interest to us. The information in the right 
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hand side array is enough to generate both v\ and v 2 and also to determine 
the type of v±. 

We can encode longer sequences of types using longer arrays, as we now 
illustrate by example: Our first gene is located at the point p = (3,4). The 
two arrays are 

II 9 23 11 14 10 

A((3,4),l,3) = 10 0; A((3,4),2,3) = -1 -1 

III 0-1-1 

This tells us that the 5 vertex types we see along the gene A are (in one of 
the two orders) 10, 14, 11, 9, 23. Looking at the first column of A((3, 4), 1, 3), 
we can see that that the vector (1,1) connects the point of type 11 to the 
point of type 9. In this way, we can draw a copy of the gene given the above 
arrays. The next routine generates the arrays for the point p = (x,y) E Z 2 . 

getOrbit (x,y ,length,epsilon) : 

1. let A be the empty array 

2. let X = x and Y — y. 3. let count= 

4. while(count<length): 

a =classif y(X,Y) 

Let M be the matrix associated to T a 

If(count= 0) then: 

append to A the (epsilon)th column of M 
if(count> 0) then: 

append to A the col. of M which does not match the last col. of A. 
Let m x and m y denote the entries of the column of M used above, 
replace X by X + m x 
replace Y by Y + m y 
increment count 

5. return(A) 

getOrbit runs until the count equals the length, and then breaks. At 
that point, the array A is returned. For our purposes, we need to get both 
sequences associated to p and k. So, we run the above routine once for e = 1 
and once for e = 2. We call the resulting pair (L, R) of arrays the k-itinerary 
of p £ Z 2 . The genes correspond to the case k — 3. 
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7.5 Checking the Dynamical Polygons 

Each array A determines a dynamical polygon P = P(A), as in §4.2. Here 
we explain how to check directly that a point s6T 2 lies in the closure of P 
or else in the interior of P (depending on our interest.) The point s is always 
given as the projection to T 2 of a certain lift s 6 R 2 which we give the same 
name. Let Pj denote the jth tile in the partition V of T 2 . As above, we 
identify Pj with a particular lift to R 2 . The following routine returns a 1 
provided that s is contained in the interior of P. 

matchOrbitOpen(s,A): 

let S — s 

loop for % = 1 to length of A 

check that isLatticeContainedOpen(S',P a J= 1. 
if false return(O) 
if true then: 

let z =getMap(a; i ,y i ) 

let S = S + z (act on S by the dynamical translation) 
endloop 

As a variant, we check that s is contained in the closure of P by using 
isLatticeContainedClosed in place of isLatticeContainedOpen. 

We can now explain exactly how we check our guess Pq for the 75 dy- 
namical polygons associated to the 75 genes. Here is the 3-step process. 

• For the jth gene Aj we let (xj,yj) be the location of the center dj of 
Ay We use get0rbit(xj,y 3 ) to generate the 3-itinerary A = (Lj,Rj). 

• We check that each Integer-Complex v representing a vertex of Pq is 
contained in the closure of P = Pa using matchOrbitClosed(f ,Lj) and 
then matchOrbitClosed(f ,Rj). 

• For each edge e of Pq we produce a IntegerComplex v contained in 
the interior of e by applying interpolate to the two endpoints of e. 
Given v we then check that one of the two runs matchOrbitOpen(f ,Lj) 
and matchOrbitOpen(t> ,Rj) returns a 0. This is to say that v is not 
contained in the interior of P. 

We list the dynamical polygons P , P 74 in §8.21 We list the centers for 
our genes in §8.31 
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7.6 The Shadowing Property 

We begin by recalling the basic setup. Let A be a gene from our inflation 
generator and let B be the gene core. Let a be the center vertex of A. Let 
A' be the strand which shadows $(A). Recall that $ is dilation by <p 3 . Let 
P be the dynamical polygon associated to A and let P' be the dynamical 
polygon associated to A'. Let (5 be the inflation map and let 7 : P -> T 2 be 
the associated special similarity. We want to verify that j(P) C P'. 

The first main task is to compute the dymamical sequences associated to 
A' We store A' as a triple {(xj,yj)} of points. Here (£2,3/2) is the point a' 
that comes close to $(a), and (xi,J/i) and (£3, 2/3) are the two endpoints of 
A'. In §8.41 we give the list of triples. 

A variant of getOrbit recovers the sequence of types from the triple. 
This variant starts at a' and moves out in either direction until it encounters 
the two endpoints. In other words, we just replace line 4 of getOrbit with a 
check that the current point equals neither (xi,yi) nor (x 3 ,y 3 ). Technically 
we could avoid this reconstruction problem just by saving the sequences 
attached to A[, ...,A' 75 , but we prefer not to store so much data. Below we 
will list out the 75 triples. 

Now we explain how to compute 7. We have already defined the routine 
psi above, which chooses some lift of \P(a) to R 2 . Here we explain an im- 
proved version, which chooses the lift of ^(a) that is contained in P. The 
idea is to find a "correction vector" A G {— 7, 7} 2 such that psi(a)+A G P. 
This always works, and we call the result canonicalPsi. It is a better 
lift of From the way we have constructed A', we know that 7 maps 
canonicalPsi(a) into P'. Thus, we can use our knowledge of canonicalPsi(a) 
to compute 7. We seek a pair of integers (m, n) such that 

T = 7o + [(m</r 4 + n0- 1 ,O)]; 7o(z, v) = h4>~ 3 x, -0" 3 y). (52) 

We set y =canonicalPsi(a) and then test small integer choices of m and n 
until we find a choice which leads to 7(1/) G P' . To test that 7(1/) G P' we 
don't need to compute P' . We simply check that j(y) follows the dynamical 
sequences associated to A'. In other words, we just go back to the definition 
of P' . We make the same kind of loop as in the routine canonicalPsi, with 
\m\, \n\ < 7. This always works. 

Once we have 7 we list out the vertices of P as P\, P\. and check that 
7 (Pi) follows the dynanical sequences associated to P'. Here k < 5 for each 
of the 75 choices. 
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7.7 Checking Coherence 

Let gi, (71024 be the first 1024 points along r . The point gj is the center of 
a gene Gj. Let j3j be the inflation map we associate to Gj using our inflation 
structure \- The list Pi, ...,/3io24 consists of multiple references to the same 
75 inflation maps. For each index j we produce a point g'j G r , always 
within 3 units of $(<?j), such that 

*M = T°%)- (53) 
This is to say that Pj(gj) = g'j. We check Equation ESI as follows: 

verif yCoherence : 
loop from j = 1 to 1024. 

1. Let fc = (the gene type of Gj) 

2. let yj =canonicalPsi(g J ) 

3. let 7 be the map computed in Equation [52] for the index k. 

4. Check that j(yj) and psi(^) agree up to a vector in {—7, 7} 2 . 
endloop 

Line 4 really says that j(yj) = ^f(gj) in T 2 . 

As a final detail, we check visually that g[, g' 1024 lie in order on r . 
This verifies the coherence of our inflation generator. 

In the interest of space, we don't list the pairs (gj,g'j) in the appendix. 
The interested reader can push the button shadow points on the Arith- 
metic Graph Control Panel of Billiard King and see these points plotted. 
The first few pairs (gj,g'j) are 

((1, 1), (3, 4)); ((1, 2), (3, 9)) ((1, 3), (3, 13)) ((2, 4), (8, 18)). 

For reference, the following routine generates g'j given gj. 

generateShadow(j) : 

1. Let k = Ij G {1, 75} be the type of the gene Gj 

2. Let a>k be the center of the kth gene in our inflation generator 

3. Let a' k be the point of A' k that shadows $(a&) 

4. Search for a point g'j G Z 2 within 3 units of &{gj) such that g'j and a' k are 
the centers of equivalent genes. 
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7.8 Checking the Partition 

Here we explain the calculations from the end of §6.61 



7.8.1 Explicit Formulas 

Let Ej be the extension of Ej to Tg guaranteed by the Extension Lemma. 
We can rigorously compute the formulas for these maps just by making a few 
calculations. (We checked the formulas on millions of points.) 

One general feature of the map Ej is that the corresponding "undefined 
set" Yj consists of two parallel 3-tori. The complementary region T 8 4 — Yj 
consists of two regions. This feature is not necessarily clear from our descrip- 
tion below, because we describe our maps in terms of coordinates on the set 
(— 4,4) 4 , which is the interior of a fundamental domain for Tg. The action 
on T 4 is obtained by piecing together the definitions across the boundaries 
of [— 4,4] 4 . In all cases, we will list the following data: 

• The linear part L of Ej. 

• A determiner function d : (— 4, 4) 4 — > R. 

• A partition of R either into 3 or 5 intervals. In the 3-interval case, the 
dividing points are and 4 and the intervals are (— oo, 0) and (0, 4) and 
(4, oo). In the 5-interval case, the dividing points are —8, —4, 0, 4 and 
the intervals are determined similarly. 

• For each interval / of the partition we give an array of the form 

V (^J^ — an fl21 fl 31 °41 
L Q-12 a 22 a 32 a 42 . 

This array stands for the vector 

V 2 ' 2 ' 2 ' 2 J 

In case all the entries of the array are we will save space by just writing [0] 
in place of the array. Our map is then given by E(x) = L(x) + v(I), where 
/ is such that d(x) G /. This definition makes sense except for the points x 
where d(x) lies in one of the dividing points listed above. 
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E 1 : L x (x) = 



1111 
0-11 



x 2 /(j) + x 3 + X A 
x 2 /<t> J 

1 1 



di(x) — xi + 1; 



[0] 



1 1 
1-1 



(54) 



R 



2 • 



L 2 (x) 



Xi 
X 2 
X 3 




2 



10 1 
1 



[0] 







d 2 (x) = x 2 + rr 3 - 



\-Xi +X 2 /(f) + X 3 <f)J 

10 1 



" 
-2 



(55) 



E% 



"0 

2 



L 3 (x) 




-2 -2 



- x 2 - x 3 + x 4 \ 

Xi/(f> — X3 + £40 

+ X 2 + 2X3 — ^40 
-Xi/(j) + X 2 + X 3 - X 4 /(j) ) 



1 

1 -1 -1 



[0] 



1 1 

-1 -1 



d 3 (x) = xi + x 4 




1 





-2 





-2 



1 _ 

02' 

0" 

2 2 
(56) 



£4 : L 4 (x) = 



"10 1 
.0 1 



/ x, \ 

2^2 

X 1 /(j) 

V rr 4 / 








[0] 



d 4 (x) = x + 1 



1 



0" 



_0 -1 0. 



(57) 
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7.8.2 Reducing to Planar Polygons 

The polygons in the partition V are subsets of the torus T|. In this section 
we explain how we arrange our computation so that we just have to consider 
polygons which are subsets of R 4 . 
Let 

T\ = R 2 /(AZ 2 ). 

Then T 4 is a square torus which is naturally a 16-fold cover of T 2 . Let V 
denote the lift of V to Trf . If P k is a polygon of V, then P k is the projection 
to T 2 of the convex hull of certain vertices {v kl: ...,v kn }. For each (ei,e 2 ) G 
{0, 1, 2, 3} 2 we can form the convex polygon P k (e 1: e 2 ) with vertex list 

{v k i + (ei, e 2 ), Ufo, + (ei, e 2 )} (58) 

and project it into T 2 . Call the resulting polygon Pfc(ei,e 2 ). The union of 
the polygons -Pfc(ei,e 2 ), taken over all indices, gives us the partition V. So, 
just to be clear, P k (e 1: e 2 ) is some concrete lift to R 2 of a polygon of the 
partition P. 

We are interested in the polygons of the partition V, which is a subset 
of the different torus T|. We now explain how to translate between V and 
V . The map 6 8 (x) = 86 (x/ 8) computes x mod 8Z. This map is the building 
block for our map ip. We introduce the maps /x_ and defined by 

H-(x, y) = ((6 8 (2x + 1), 6 8 (2x - 1), 6 8 (2y + 5 8 (2y - 

fi + (x, y) = ((6 8 (2x - 1), 6 8 (2x + 1), 5 8 (2y - 5 8 (2y + 0" 1 )) 
We also introduce the map 

$(x) = (6 4 (x/(2<f))),6 4 (x/2)) = (46(x/(8cf ) )),46(x/8)) e T 2 . (59) 

We compute that 

/!+ o ip(x, 1) = ip(x, 1); /i_ o —1) = ip(x, —1). 

It follows from this last equation that ji + maps the polygons of V to the 
polygons of V which partition T 2 (+). A similar statement holds for //_. 

We can think of /i + and as being defined on i? 2 , and thus the polygons 
of V all have the form 

At±(Pfe(ei,e 2 )) 
for fc = 1, ...,26 and G {0, 1, 2, 3}. 
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7.8.3 The Calculation 



Here's the routine that produces a lift of the 02 nd point of fi±(Q), where 
Q — Q ^ (0,3, £14). The integer 05 G {0, 1} toggles /i+ and 

getVertexPlanar(a! ,02,03,04,05): 

Let v = x + iy be the a 2 vertex of our particular lift of P ai . 
Let X2 = x\ + a 3 and yi — y\ + 04. (Equation [581 ) 

If 05 = then return (2xi + 1, 2x2 — 1, 2?/2 + 1/0, %2 — V0)- ( a PPly /•*+•) 
If 05 = 1 then return (2xi — 1,2x2 + 1,27/2 — 1/0,2^2 + 1/0)- (apply A 4 --) 

Next, we have a routine dec8, which computes x mod 8Z, where x G 
Here dec8 works just like our routine dec above, and is checked in the same 
way. In brief dec8(x)=8dec(x/8). 

Our routine subdivide starts with a polygon P C R 4 and adds one 
new point interpolate^ ,t>j + i) between consecutive vertices Vi and v i+ i of 
P. This, the new polygon has twice as many vertices as the old one. Our 
routine getPolygon subdivides a polygon in R 4 ten times and then projects 
it into Tg using dec8 componentwise on all the vertices. 

We use a routine getTracePoint to get an interior point of Q. This rou- 
tine returns the point interpolate(x 3 , interpolate^!, X2)) where Xi, X2, X3 
are the first three vertices of Q. 

For each of the relevant indices (ai, a 3 , a 4 , a 5 ) we let P be the poly- 
gon produced by getPolygon and we let x G P be the point produced 
by getTracePoint. Our routine getltinerary computes the itinerary of 
x in a straightforward way, just by computing the orbit and checking which 
component Ck{^k) contains the relevant point x& of the orbit for k — 1, ...,8. 
The itinerary is then e±, e§, as in §6.61 Finally, for each vertex y of P we 
perform the following routine 

verif yltinerary(e,?/): 
Let yi = y. loop from k = 1 to 8. 
Verify that y k G C k (e k ). 

If false then return(f alse). Otherwise continue. 
Let y k +i = Ekiyk) using the e k extension of E k . 
return(true) 
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8 Appendix 

8.1 The Polygons in the Partition 

Here we list out the coordinates of the polygons in our partition V of T 2 . 
These polygons are drawn in Figure 2.5. Each vertex of Pj has the form 



fa + ai<f> a 2 + a 3 0\ 



We will simply list such a vertex as a row 

a ai a 2 a 3 

in an array whose other rows correspond to the other vertices of P, 



3 ' 



PI P2 
5 -3 -5 3 



P3 





10 -6 -6 4 -10 6 6 -4 5 3 

2 -1 -2 1 ^ ~1 'I \ 

o —6 — 6 4 



PS P6 

0000 000 

13 -8 -6 4 -13 8 6 -4 

5 -3 -6 4 -5 3 6 -4 



P9 P10 

3 -2 -4 2 -3 2 4 -2 

-5 3 2 -2 5 -3 -2 2 

3 -2 -2 -3 2 2 

-2 1 -2 2 -1 2 

6 -4 -8 4 -6 4 8 -4 



P12 P13 P14 

10 -6 -6 4 -10 6 6 -4 

-5 3 4 -2 5 -3 -6 4 -5 3 6 -4 

3 -2 -3 2 4 -2 3 -2 -4 2 



















-5 


3 










-10 


6 


6 


-4 




-5 


3 


6 


-4 


P8 










-9 


5 


4 




-4 


-4 


2 


-2 







-1 





-2 







-14 


8 


4 




-4 



P7 

9 -5 -4 4 
4-220 

10 2 
14 -8 -4 4 



P16 

-13 8 6 -4 

-5 3 6 -4 

3 -2-4 2 

-5 3 2 -2 



P17 P18 

12 -8 -8 4 -12 8 8 -4 

7 -5 -8 4 -7 5 8 -4 

4 -3 -4 2 -4 3 4 -2 



P20 P21 P22 

4 -2 2 1 -1 -1 1 

-1 1 2 6 -4 -6 4 -6 4 6 -4 

-4 3 6 -2 -2 1 2 -1 















5 


-3 


-4 


2 


-3 


2 








P15 








13 


-8 


-6 


4 


5 


-3 


-6 


4 


-3 


2 


4 


-2 


5 


-3 


-2 


2 


P19 








-4 


2 


-2 





1 


-1 


-2 





4 


-3 


-6 


2 


P23 








6 


-3 


-2 


2 


-2 


2 


2 





-4 


3 


4 


-2 


4 


-2 









1 


-1 








-1 


1 








P26 








-2 


1 








2 


-1 




















3 


-2 


-4 


2 


-3 


2 


4 


-2 


5 


-3 


-4 


2 


1 


-1 


-2 





-1 


1 


2 











-2 





-1 





-2 





1 





2 





-5 


3 


2 


-2 


4 


-3 


-4 


2 


-4 


3 


4 


-2 
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8.2 The Dynamical Polygons 

Here we list the 75 dynamical polygons associated to the 75 genes. We use 
the same notation as above. 



PC) 

13 -8 

-21 13 16 -10 

13 -8 -10 6 
5-3 



PI 

-15 9 8 -6 

19 -12 -8 4 

11 -7 -8 4 

19 -12 -18 10 



P2 

9 -5 -6 4 

17 -10 -6 4 

-17 11 10 -6 

17 -10 -16 10 



P3 

13 -8 -6 4 

5-3-6 4 

13 -8 -16 10 

-21 13 10 -6 



P4 








16 


-10 


-8 


4 


3 


-2 


-2 





-10 


6 


-2 





24 


-15 


-18 


10 


P8 








-10 


6 


12 


-8 


3 


-2 


-4 


2 


16 


-10 


-14 


8 


24 


-15 


-14 


8 



P5 

14 -8 -4 4 

1 2 

-12 8 2 

22 -13 -14 10 



P9 

-20 13 14 -8 

1 0-2 2 

-12 8 14 -8 



P6 








19 


-11 


-8 


6 


-15 


10 


18 


-10 


-7 


5 


8 


-4 


-15 


10 


8 


-4 



P10 

-18 11 6 -4 

3 -2 -10 6 

-10 6 6 -4 



P7 

10 -6 -6 4 

5 -3-6 4 

-3 2 4 -2 



Pll 

23 -14 -16 10 

10 -6 
2-100 



P12 

-16 10 14 -! 

5 -3-2 2 

26 -16 -12 8 

5 -3 -12 8 



P13 

-17 10 4 -4 

-4 2 4 -4 

17 -11 -12 6 



P14 








-12 


8 


14 


-8 


-20 


13 


14 


—8 


14 


— 8 


-12 


8 


22 


-13 


-12 


8 



-15 


10 


18 


-10 


19 


-11 


-8 


6 


-2 


2 


2 






P16 

-23 14 16 -10 

-10 6 

11 -7 -10 6 



P20 

24 -15 -18 10 

3 -2-8 4 

16 -10 -8 4 



P17 








16 


-10 


-14 


8 


-5 


3 


2 


-2 


-26 


16 


12 


-8 


-5 


3 


12 


-8 


P21 








1 





-4 


4 


22 


-13 


-14 


10 


14 


-8 


-4 


4 



P18 



17 


-10 


-4 4 


4 


-2 


-4 4 


-4 


3 


6 


P22 






-19 


12 


8 -4 


-6 


4 


8 -4 


2 


-1 


-2 2 



P19 

13 -8 

-8 5 10 -6 





P23 

9 -5-2 2 

-12 8 14 -8 

-4 3 4 -2 



P24 








-5 


3 


4 


-2 


16 


-10 


-6 


4 


-5 


3 


-6 


4 


-26 


16 


20 


-12 



P25 

19 -12 -16 10 
6-400 

-15 9 10 -6 



P26 P27 

17 -10 -16 10 13 -8 -16 10 

-17 11 10 -6 -21 13 10 -6 

4 -2 



P28 

13 -8 

-21 13 16 -10 





P29 

19 -12 -8 4 

6 -4 -8 4 

-15 9 8 -6 



P30 

-9 5 2 -2 

12 -8 -14 8 

4-3-4 2 



P31 








5 


-3 


-4 


2 


-16 


10 


6 


-4 




-3 


6 


-4 


26 


-16 


-20 


12 



P32 








P33 








P34 






P35 






-6 


4 








4 


-3 


-2 





24 


-15 


-14 8 


14 


-8 


-12 8 


-19 


12 


16 


-10 


-17 


10 


14 


-10 


16 


-10 


-14 8 


22 


-13 


-12 8 


2 


-1 








-4 


2 


-2 





3 


-2 


2 -2 


1 





4 -2 



P36 

-17 10 4 -4 

-4 2 4 -4 

4 -3-6 2 



P37 

-13 8 

8 -5 -10 6 





P38 








-16 


10 


8 


-4 


-3 


2 


2 





10 


-6 


2 





-24 


15 


18 


-10 



P39 








-14 


8 


4 


-4 


-1 





-2 





12 


-8 


-2 





-22 


13 


14 


-10 



P40 






P41 






P42 






P43 








-10 


6 


6 -4 


-24 


15 


14 -8 


26 


-16 


-20 12 


-19 


12 


16 


-10 


-5 


3 


6 -4 


-11 


7 


14 -8 


5 


-3 


6 -4 


-6 


4 








3 


-2 


-4 2 


10 


-6 


-12 8 


-8 


5 


6 -4 


15 


-9 


-10 


6 
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P44 








P45 






P46 








P47 








-13 


8 


16 


-10 


-24 


15 18 


-10 


-1 





4 


-4 


19 


-12 


-8 


4 


21 


-13 


-10 


6 


-3 


2 8 


-4 


-22 


13 


14 


-10 


6 


-4 


-8 


4 














-16 


10 8 


-4 


-14 


8 


4 


-4 


-2 


1 


2 


-2 



P48 








P49 








P50 








P51 






18 


-11 


-10 


6 


-10 


6 


-2 





-12 


8 


2 





24 


-15 


-14 8 


5 


-3 


6 


-4 


24 


-15 


-18 


10 


22 


-13 


-14 


10 


11 


-7 


-14 8 


-16 


10 


6 


-4 


11 


-7 


-2 





9 


-5 


2 





-10 


6 


12 



P52 






P53 








P54 








P55 








9 


-5 


-12 8 


-26 


16 


20 


-12 


6 


-4 








-4 


3 


2 





22 


-13 


-12 8 


-5 


3 


-6 


4 


19 


-12 


-16 


10 


17 


-10 


-14 


10 


-12 


8 


14 -8 


8 


-5 


-6 


4 


-2 


1 








4 


-2 


2 






P56 



-24 


15 


14 


-8 


-16 


10 


14 


-8 


-3 


2 


-2 


2 


P60 








10 


-6 


2 





-24 


15 


18 


-10 


-11 


7 


2 






P57 

-13 8 

21 -13 -16 



P61 

12 -8 -2 

-22 13 14 -10 

-9 5 -2 



P58 



-19 


12 


8 -4 


-6 


4 


8 -4 


15 


-9 


-8 6 


P62 






10 


-6 


-12 8 


-3 


2 


4 -2 


-16 


10 


14 -8 


-24 


15 


14 -8 



P59 








-18 


11 


10 


-6 


-5 


3 


-6 


4 


16 


-10 


-6 


4 


P63 








18 


-11 


-6 


4 


-3 


2 


10 


-f 


10 


-6 


-6 


4 



P64 








-23 


14 


16 


-10 


-10 


6 








-2 


1 









P65 






P66 






P67 






-13 8 








15 


-9 -8 


6 


-13 8 


6 


-4 


21 -13 


-16 


10 


-19 


12 8 


-4 


-5 3 


6 


-4 


-13 8 


10 


-6 


-11 


7 8 


-4 


-13 8 


16 


-10 


-5 3 








-19 


12 18 


-10 


21 -13 


-10 


6 



P68 








P69 








P70 








P71 






-26 


16 


12 


-8 


17 


-10 


-4 


4 


— 5 


3 


12 


-8 


14 


-8 


-12 8 


-5 


3 


12 


-8 


4 


-2 


-4 


4 


-18 


11 


12 


-8 


-20 


13 


14 


8 


-5 


-4 


2 


-17 


11 


12 


-6 


16 


-10 


-14 


8 


-7 


5 


14 -1 






P72 








P73 








P74 














26 


-16 


-12 


8 


23 


-14 


-16 


10 


5 


-3 


-12 


8 








5 


-3 


-12 


8 


10 


-6 








18 


-11 


-12 


8 








-8 


5 


4 


-2 


-11 


7 


10 


-6 


-16 


10 


14 


-8 





8.3 The Gene Locations 

Here we list the coordinates for the points ao, ...,074. The point aj is 
central vertex of the gene Aj. 



aO 


al 


a2 


a3 


a4 


a5 


a6 


al 


a8 


a9 


3 


4 


4 


4 


4 


4 


3 


5 


6 


6 


4 


5 


6 


7 


10 


11 


13 


15 


16 


17 



alO 


all 


al2 


al3 


al4 


al5 


al6 


al7 


al8 


al9 


7 


8 


10 


11 


14 


16 


16 


14 


13 


12 


19 


16 


17 


16 


20 


31 


33 


32 


33 


31 



62 



a20 


a21 


a22 


a23 


a24 


al5 


a26 


a27 


a28 


a29 


17 


17 


15 


14 


13 


12 


12 


12 


16 


17 


49 


50 


50 


49 


50 


50 


51 


52 


64 


65 


a30 


a31 


a32 


a33 


a34 


a35 


a36 


a37 


a38 


a39 


18 


19 


20 


20 


22 


22 


32 


33 


33 


33 


66 

\J\J 


65 


65 


64 


64 


65 


79 


81 


78 


77 


a40 


a41 


a42 


a43 


a44 


a45 


a46 


a47 


a48 


a49 


32 


31 


32 


41 


41 


41 


41 


43 


45 


46 


73 


72 


70 


73 


71 


68 


67 


67 


67 


68 


a50 


a51 


a52 


a53 


a54 


a55 


a56 


a57 


a58 


a59 


46 


69 


69 


68 


67 


67 


65 


63 


62 


55 


69 


137 


138 


139 


139 


140 


140 


137 


136 


142 


a60 


a61 


a62 


a63 


a64 


a65 


a66 


a67 


a68 


a69 


54 


54 


52 


51 


50 


63 


62 


62 


61 


55 


141 


140 


135 


132 


135 


213 


212 


210 


207 


214 



a70 


all 


a72 


a73 


a74 


48 


48 


73 


76 


86 


210 


211 


282 


275 


279 



8.4 The Inflation Data 

Here we list the data for the strands A' Q , ...,A' 74 that are associated to the 
genes A , A 7i . The listing 

A, 
x 

yo 

Xi 

Vi 

x 2 

V2 



63 



indicates that the two endpoints of Aj are (x , y ) and (x 2 , y 2 ), and the point 
(xi, y±) is the one lying within 3 units from the center point of D(Aj). (There 
might be several such points, but we make some choice in each case.) Here 
is the listing. 





A 




Ai 


A 2 


A3 


A 4 


A 5 


A 6 


A 7 


A 8 




A g 


A10 


An 


A12 






9 




11 


15 


17 


12 


15 


12 


15 


21 




25 


24 


33 


37 






17 




16 


22 


25 


39 


43 


50 


64 


63 




68 


76 


73 


72 






11 




16 


17 


16 


16 


17 


12 


21 


25 




24 


31 


34 


41 






17 




22 


25 


30 


43 


46 


56 


63 


69 




72 


80 


68 


72 






15 




17 


16 


12 


17 


12 


12 


25 


24 




24 


34 


37 


41 






22 




25 


30 


34 


46 


50 


60 


68 


71 




76 


76 


72 


67 




A13 




A 14 




A15 


Ai6 


A17 


Ais 


A19 


A20 




A21 




A22 


A23 


A24 


A25 


41 




59 




67 


67 


62 


58 


54 


67 




71 




66 


63 


58 


51 


67 




81 




128 


136 


136 


140 


136 


204 




208 




217 


208 


212 


208 


46 




58 




67 


67 


59 


54 


51 


71 




71 




62 


58 


54 


50 


69 




85 




132 


140 


137 


140 


132 


208 




213 




211 


208 


212 


211 


45 




58 




67 


62 


58 


54 


46 


71 




66 




63 


58 


51 


51 


71 




89 




136 


135 


140 


136 


136 


212 




217 




208 


212 


208 


216 


A 2 6 




A 2 7 




A 2 8 


A 2 9 


A30 


A31 


A32 


A33 




A34 




A35 


A 36 


A37 


A 38 


52 




51 




64 


67 


72 


78 


85 


85 




89 




92 


130 


134 


143 


211 




216 




271 


272 


279 


275 


279 


276 




267 




271 


334 


338 


335 


51 




50 




68 


71 


75 


79 


84 


86 




92 




92 


135 


138 


139 


216 




221 




272 


276 


280 


276 


276 


270 




271 




276 


336 


344 


331 


50 




46 




71 


72 


76 


84 


84 


89 




92 




92 


134 


143 


139 


221 




225 




276 


279 


275 


279 


271 


267 




275 




280 


338 


340 


326 


A 3 9 




A40 




A41 


A42 


A43 


A44 


A45 


A 4 6 




A47 




A 4 8 


A49 


A50 


A51 


140 




140 




135 


132 


174 


173 


177 


172 




178 




187 


189 


193 


288 


331 




310 




309 


297 


313 


305 


293 


288 




280 




284 


283 


289 


576 


141 




135 




130 


134 


173 


175 


172 


175 




181 




189 


194 


195 


291 


325 




310 




305 


297 


310 


300 


287 


283 




284 




284 


289 


292 


580 


143 




130 




130 


139 


173 


177 


173 


178 




181 




193 


195 


190 


291 


322 




305 




301 


300 


305 


297 


284 


280 




288 




289 


292 


296 


584 


A52 




A53 




A54 


A55 


A 5 6 


A57 


A 5 8 


A59 




A 6 o 




A 6 l 


A 6 2 


A 6 3 


A 64 


291 




291 




284 


285 


279 


271 


266 


237 




232 




229 


224 


219 


211 


580 




589 




585 


588 


597 


580 


581 


600 




602 




598 


576 


563 


568 


291 




287 




283 


284 


274 


265 


261 


231 




228 




230 


219 


216 


211 


585 




589 




588 


593 


593 


581 


575 


602 




598 




592 


572 


559 


572 


291 




284 




284 


279 


274 


261 


262 


228 




228 




232 


219 


211 


206 


589 




585 




593 


597 


589 


575 


572 


598 




593 




589 


568 


563 


567 



A 6 5 


A 66 


A 6 7 


A 6 8 


A 6 9 


A70 


A71 


A72 


A 73 


A 74 


271 


266 


262 


261 


236 


206 


203 


303 


321 


359 


902 


903 


894 


877 


907 


890 


890 


1195 


1170 


1182 


265 


261 


264 


258 


233 


203 


202 


308 


322 


363 


903 


897 


889 


878 


906 


891 


894 


1195 


1165 


1182 


261 


262 


267 


257 


231 


202 


202 


308 


325 


363 


897 


894 


885 


881 


903 


894 


898 


1190 


1169 


1177 



64 



9 References 

[B] P. Boyland, Dual Billiards, twist maps, and impact oscillators, Nonlin- 
earity 9 (1996) 1411-1438 

[D], R. Douady, These de 3-eme cycle, Universite de Paris 7, 1982 

[DT] F. Dogru and S. Tabachnikov, Dual Billiards, Math Intelligencer vol. 
27 No. 4 (2005) 18-25 

[G] D. Genin, Regular and Chaotic Dynamics of Outer Billiards, Penn State 
Ph.D. thesis (2005) 

[GS] E. Gutkin and N. Simanyi, Dual polygonal billiard and necklace dy- 
namics, Comm. Math. Phys. 143 (1991) 431-450 

[Ke] R. Kenyon, Inflationary tilings with a similarity structure, Comment. 
Math. Helv. 69 (1994) 169-198 

[Ko] Kolodziej, The antibilliard outside a polygon, Bull. Polish Acad Sci. 
Math. 37 (1989) 163-168 

[M] J. Moser, Stable and Random Motions in Dynamical Systems, with Spe- 
cial Emphasis on Celestial Mechanics, Annals of Math Studies 77, Princeton 
University Press (1973) 

[N] B.H. Neumann, Sharing Ham and Eggs, 

summary of a Manchester Mathematics Colloquium, 25 Jan 1959 

published in Iota, the Manchester University Mathematics students' journal 

[T] S. Tabachnikov, Geometry and Billiards, A. M.S. Mathematics Advanced 
Study Semesters (2005) 

[VS] F. Vivaldi, A. Shaidenko, Global stability of a class of discontinuous 
dual billiards, Comm. Math. Phys. 110 (1987) 625-640 



65 



